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ABSTRACT 


The identification and control of vapor phase reaction 
kinetics to produce pigments by homogeneous nucleation have 
been achieved. A vapor phase apparatus has been designed, 
fabricated, and calibrated through l800°C. Vapor phase re- 
actions have been analyzed, calculations have been made, 
and powders of alumina, rutile, zinc orthotitanate (in a 
mixed phase) calcium tungstate, and lanthana have been pro- 
duced by homogeneous nucleation. Electron microscopy shows 
uniform particle morphology and size and supports anticipated 
advantages of vapor-phase homogeneous nucleation; namely, 
purity, freedom from defects, and uniform particle sizing 
without grinding. 


descriptors: 

Thermal control space -stable pigments. Vapor phase, 
Nucleation, Alumina, Rutile, Zinc Orthotitanate, 
Lanthana, Calcium Tungstate, Ceramics, Gas Mixing 


ii 



FOREWORD 


This report was prepared by The Ohio State University Research 
Foundation, under MSA Contract No. NAS8-21317, entitled, "Prepara- 
tion of Pigments for Space-Stable Thermal Control Coatings," This 
reports covers the period from June 1, 1968 through March 31 > 1972. 

Dr. William B. Campbell was the principal investigator for this 
program. He wtoS ably assisted by R. G. Smith, J. K. Cochran, J. W. 
Hinton, S. G. Nychas , J. E. Seaver, J. W. Randall, R. J. Versic 
and J. E. Burroughs. Full or partial support from this program 
provided the means for R. G. Smith, J. K. Cochran, J. W. Hinton, 

S. G. Nychas and R. J. Versic to fulfill the requirements for the 
Doctor of Philosophy degree. The Master of Science degree was re- 
ceived by J. W. Randall. 

This work was performed under the Space Sciences Laboratory of 
the George C. Marshall Space Flight Center. Mr. Daniel W. Gates 
was the Contracting Officer's Technical Representative and the 
Project Manager. 


Respectfully submitted 
OSU Research Foundation 


William B. Campbell, Ph.D. 
Associate Professor 


Approved: 



OSU Research Foundation 


iii 



SUMMARY 


The development of vapor phase technology in the preparation of 
space-stable pigments has been demonstrated by the production of 
very fine powders with des crete morphology and an average particle 
size in the range of 0.35 to. 2.5 p. 

The critical parameters for the controlled production of rutile 
and alumina have been used as the basis for the development of mono- 
sized particles of zinc orthotitanate. Zinc orthotitanate (Zn2Ti0 4 ) 
has been found to be one of the better space-stable pigments for use 
in thermal -control paints of low o: and low Q!/e ratio (solar absorpt- 
ance and infrared emittance). Because of these stringent require- 
ments, the development of vapor phase reactions to produce calcium 
tungstate and lanthana by homogeneous nucleation was included in 
this program. 

Vapor phase apparatus has been designed, fabricated, and cali- 
brated. Vapor phase reactions have been analyzed, calculations 
have been made, and powders of alumina, rutile, zinc orthotitanate 
(in a mixed phase), calcium tungstate and lanthana have been pro- 
duced by homogeneous nucleation. Submicron size particles required 
for efficient scattering of solar radiation in the 0.35 to 2.5 P 
wavelength range have been produced. Excellent sub-angular mor- 
phology was exhibited in the electron photomicrographs. 
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PREPARATION OF PIGMENTS FOR SPACE-STABLE THERMAL CONTROL COATINGS 


I . INTRODUCTION 

During the past decade, efforts have been directed toward producing 
thermal control coatings which are stable to solar radiation in space. 
Inorganic pigments, once believed to be extremely stable, have been 
shown to be subject to ultraviolet and proton degradation. The degra- 
dation mechanism appears related to the strong absorption of high-energy 
ultraviolet radiation and may result from a defect structure in the 
pigment crystals. Crystal size, structural integrity, and surface con- 
dition apparently influence the performance and degradation resistance 
of a pigment. 

The primary purpose of a thermal control coating is to reflect as 
much radiation as possible. Maximum reflectance is obtained when the 
particle size of the pigment is of the same order of magnitude as the 
wavelength of the incident radiation. For solar radiation, the maximum 
intensity is between 0.35 and 2.5 microns. 

In the past eight years, the major emphasis in pigment research 
for thermal control coatings has been on zinc oxide which has consist- 
ently shown the greatest stability to ultraviolet irradiation in a 
vacuum. Of the new pigments under investigation, zinc orthotitanate 
has shown considerable promise as a new, potentially stable white 
pigment. 

The first zinc orthotitanate (Zn 2 TiCu) pigment powder was prepared 
by mixing ZnO and Ti0 2 (anatase) in the mole ratio (2:1). Residual ZnO 
was removed by washing with acetic acid (complete removal is seldom 
achieved) . Stoichiometric considerations suggest that residual Ti0 2 
remains, that other titanates are formed, or that titanium is present 
in interstitial solid solution in the orthotitanate lattice. 1 Subse- 
quent, solid state reactions have reportedly reduced the stoichiometric 
effects. 2 

Zinc orthotitanate is not completely space-stable and will degrade 
under space simulation. One of the primary factors believed responsible 
for promoting damage to the pigment powder is grinding; however, the 
extraction with acetic acid and washing are other important variables. 

A grinding or mulling operation is usually employed to obtain submicron- 
size particles that produce the efficient scattering of solar radiation 
in the 0.35 to 2.5p wavelength range. 

The development of vapor -phase technology in the preparation of 
space-stable pigments has been demonstrated by the production of rutile 
powders having discrete morphology. Powders produced by homogeneous 
nucleation in the vapor phase have controlled chemical composition; are 
monosized crystalline particles of controlled particle size; and retain 
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their structural integrity. The ability of the vapor phase process to 
control the parameters that influence reflectance provides many advant- 
ages in comparison to other powder production techniques. Systematic 
variations of particle size have been controlled with hot zone temper- 
ature, residence time, and injector design. The current investigation 
was designed to study the production of stable pigments for thermal 
control coatings. Materials include zinc titanates, calcium tungstate 
and lanthanum oxide. 

Two major limitations retarding the production of multiple oxides 
and rare earth compounds ; gas mixing and halide generation have been 
studied. A detailed mixing study has been performed and a reactor de- 
signed from the results of this study. Halide generation of multiple 
chlorides has been achieved by the direct chlorination of the multiple 
oxide. Calcium tungstate has been produced as a major phase in the jet 
mixing reactor. Zinc orthotitanate and lanthanum oxide were obtained 
as minor phases only. Reaction kinetics were limiting in production of 
multiple oxide phases. 


A. PROGRAM OBJECTIVE 

Since i960 several exploratory investigations have established the 
feasibility of production of highly reproducible materials by vapor- 
phase reactions. Also, the morphology of vapor-phase productions has 
been controlled in these rapid reaction systems. Product morphology 
has ranged from whiskers to bulk crystals to monodisperse powders. 

The objective of this investigation is to apply the vapor -phase 
technology developed in earlier investigations to the products of 
space-stable pigments of zine titanates, zirconates, etc. The primary 
material system to be considered is zinc orthotitanate; however, rutile, 
aluminum oxide, calcium tungstate and lanthanum oxide are included. 

Zinc orthotitanate (ZngTiCU), which was found to be one of the 
better space-stable pigments of low a and low a/e ratios (solar absorp- 
tance and infrared emittance), was selected as the primary material for 
production by homogeneous nucleation in a quantity sufficient for pig- 
ment evaluation. 

The vapor reaction system has the potentials of close control of 
chemical composition, production of monosized cryst all ine particles, 
and retention of structural integrity -- all with sizeable production 
batches . 
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B. SURVEY OF LITERATURE 

1. Pigment Properties for Thermal 
Control of Space Vehicles 

Temperature requirements for space vehicles impose stringent 
restrictions on the properties of thermal control coatings. Satellites 
containing electronic equipment must maintain temperatures of -10°C to 
100°C and should be kept in the temperature range of 20° to 40°C for 
reliable operation. Manned space vehicles must maintain an even more 
narrow temperature range and must not exceed 45°C for more than a few 
minutes . 3 


The average equilibrium surface temperature of a vehicle in 
space is given by the following general equation: 3 


T 


(a/e) (Ps+Pa) + Pe 


Act 


( 1 ) 


where 


a 

€ 


Ps 

Pa 

Pe 

A 

cr 


solar absorptivity of external body surface; 
hemispherical emittance of external body surface, 
assumed to be equal to the absorptance of the 
Earth-emitted radiation; 

direct solar radiation incident upon the body; 
Earth-reflected solar radiation upon the body; 
Earth-emitted radiation incident upon the body; 
surface area of the body; and 
Boltzman's constant. 


To maintain low temperatures in a space vehicle, a coating of 
low (a/e) ratio and a low a is desired, as predicted by Eq. (l). A 
number of inorganic ceramic pigments, including zinc orthotitanate, 
possess low a/e ratios. 

In recent years inorganic pigments, once believed to be 
extremely stable, have been shown to be subject to ultraviolent and 
proton degradation. Considering the mechanisms involved, pigments having 
a minimum of structural defects would be the most space stable. When 
pigments degrade, the solar absorption (as) increases, thereby in- 
creasing the a/e ratio. 

The primary purpose of a thermal control coating is to reflect 
as much radiation as possible; thereby preventing the radiation from 
reaching the interior of the space vehicle. The most important factor 
to consider in the reflectivity of a two-phase coating are (l) the 
particle size of the pigment, (2) the relative index of refraction of 
the pigment and the binder, and (3) the volume of the pigment present 
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in the coating. 4 Maximum backscattering or reflectance is obtained 
when the particle size of the pigment is of the same order of magnitude 
as the wavelength of the incident radiation. For solar radiation, the 
maximum intensity is between 0.35 and 2.5p. 

For greatest reflectance, the index of refraction of the pig- 
ment should be significantly higher than that of the binder and the 
volume of particles should be as great as possible. A pigment suitable 
for a thermal control coating should (l) have a low ( a/e ) ratio, (2) be 
stable when subjected to ultraviolent radiation in a vacuum, (3) have 
a particle size in the 0.35 to 2 . 5 m- range, and (4) have a high index of 
refraction. 


2. Phase Equilibrium in the Zn0-Ti0 2 System 

Compounds reported for the Zn0-Ti0 2 system are zinc orthoti- 
tanate (Zn 2 TiC> 4 ), zinc metatitanate (ZnTi0 3 ), and zinc sesquititanate 
(Zn 2 Ti 3 0 8 ). The study of phase equilibria in the Zn0-Ti0 2 system made 
by Dulin and Rase 5 shows decomposition of zinc metatitanate into zinc 
orthotitanate and rutile at approximately 900°C (see Fig. l) by the 
following reaction: 


2 ZnTi0 3 — * Zn 2 Ti0.j + Ti0 2 . (2) 

In the composition region Zn0-Ti0 2 to Zn0-2Ti0 2 , Bartram and 
Slepetys 6 report the formation of Zn 2 Ti 3 0 8 as the predominant phase 
from 600 to 900°C . Like the metatitanate, the sesquititanate becomes 
•unstable in the temperature region 900° - 1000°C . At temperatures 
greater than 1000°C, the sesquititanate decomposes to orthotitanate 
and rutile by the following reaction: 

Zn 2 Ti;=P 8 — » Zn 2 Ti 04 + 2Ti0 2 . (3) 

The presence of Zn 2 Ti 3 0 a may have been overlooked because of similarity 
between Zn 2 TiC >4 a*id Zn 2 Ti 3 0 a structures. Variations in X-ray data have 
been previously attributed to solid solution of Ti0 2 in Zn 2 Ti 04 instead 
of Zn 2 Ti^D 8 . One point on which there is general agreement is that zinc 
orthotitanate is the only stable zinc oxide-titania compound above 
1000°C. 


The structures of zinc orthotitanate, zinc metatitanate, and 
zinc sesquititanate are inverse spinel, hexagonal, and defect spinel, 
respectively. Zinc orthotitanate with an inverse spinel structure has 
32 oxygen atoms in an approximately cubic close-packing arrangement, 
with 8 Zn' + ions in A sites with tetrahedral coordination, and 8 Zn 2+ + 

8 Ti 4 ions distributed randomly in B sites with octahedral coordination 
to the oxygen atoms. Zinc sesquititanate crystallizes in a defect 
spinel structure with 8 Zn 2+ ions in the tetrahedral positions and 
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12 Ti 4+ ions occupying the octahedral positions. Zinc metatitanate has 
an ilmenite-type structure. 


3. Vapor Phase Powder Technology Developments 

Aluminum oxide powder was produced by Campbell 7 from homogene- 
ous nucleation in the vapor phase. Aluminum chloride, carbon dioxide, 
and hydrogen were reacted in an isothermal hot zone. Variables such as 
temperature, gas velocity, system pressure, and gas composition were 
shown to be controlling parameters for powder and whisker growth. 

Schaffer and Jones 0 produced other oxide powders (i.e., Si0 2 , 
ZrO s , Ti0 2 , and ZnO and the double oxide ZnTi0 3 ) using a halide hydroly- 
sis reaction and subsequent homogeneous nucleation of the powder in the 
gas phase using equipment similar to that of Campbell. Both chlorination 
and vaporization of zinc were used for production of ZnO and ZnTi0 3 . 
Chlorination of zinc, proved unsatisfactory because of its high vapor 
pressure. Vaporization of zinc proved difficult to control because 
other reacting gases passed over its surface and oxide formation limited 
evaporation . 

Vapor-grown zinc oxide pigmenting powders are produced in 
large quantities commercially by reaction of zinc oxide with coal at 
elevated temperatures. The coal reduces zinc oxide to zinc metal vapor, 
and the zinc vapor is later reoxidized to form zinc oxide powder by 
homogeneous nucleation. 
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II. PRINCIPLES OF VAPOR PHASE POWDER PRODUCTION 
The fundamental considerations for a vapor phase reaction system 

are 

(1) the thermodynamic behavior of the system at selected 
temperatures; 

(2) the kinetics of nucleation, reaction, and evaporation; 
and 

(3) a basis for process control. 

Formation of desired products may be accomplished by introducing 
suitable metal vapor species into a heated reaction zone at constant 
temperature and pressure followed by H 2 and C0 2 to provide the remaining 
reactants necessary to cause homogeneous nucleation. By carrying out 
the chemical reaction in a zone of constant temperature and pressure, 
the equilibrium conditions of the reaction can be predicted from thermo- 
dynamics . Appendix A, B and C contain the data, and computer source 
listings used in the following calculations. 

A. THERMODYNAMIC S OF ALUMINUM OXIDE (A1 2 0 3 ) REACTION 

Aluminum oxide powders can be prepared by the following reaction: 

2 A1C1 3 + 3 H 2 + 3 C0 2 ^ AlaOa + 6 HC1 + 3 CO (4) 

The reaction free energy and equilibrium constant over the tempera- 
ture range 500-2000°K, are summarized in Table I. 

B. THERMODYNAMICS OF TITANIUM DIOXIDE (Ti0 2 ) REACTION 

Titanium dioxide (rutile) powder has been formed by homogeneous 
nucleation according to the reaction: 

TiCU g + 2 H 2g + 2 C0 2g - Ti0 2g + 4 HCl g + 2 C0 g . (5) 

The free energy change over the temperature range of 300 to 2000 °K 
reported in Table II. 
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Table I 

- Free -Energy Changes 
Reactions 

and Equilibrium Constants 

for Aluminum Oxide 

.T 

(°c) 

T 

(°K) 

AG° 

Reaction* 

(kcal/mole) 

% 

Reaction* 

227 

500 

-63.3 

3.3 x 10 27 

727 

1000 

-75-8 

3.0 x 10 16 

927 

1200 

- 80.2 

3.4 x 10 14 

1127 

1400 

-85.4 

1.8 x 10 13 

1227 

1500 

-87.6 

5.0 x 10 12 

1327 

1600 

-90.5 

2.0 x 10 12 

1427 

1700 

- 92.6 

6.9 x 10 11 

1527 

1800 

-94.6 

2.7 x 10 11 

1627 

1900 

-97.4 

1.4 x 10 11 

1727 

2000 

-100.0 

7.4 x IO 10 


*2 A1C1 3 + 3H 2 + 3 CO 


A1 2 0 3 + 6 HC1 + 3 CO 



Table II - Free-Energy Changes and Equilibrium Constants as 
a Function of Temperature for Titanium Dioxide 
Production 


T 

(°c) 

T 

(°k) 

AG° 

Reaction* 

(kcal/mole) 

Reaction* 

27 

300 

- 4.60 

2.25 x 10 3 

127 

400 

-10.55 

5.82 x 1CP 

227 

500 

-17.15 

3.14 x id 7 

327 

600 

. -18.20 

4.27 x 1CP 

427 

700 

-21.85 

6.64 x 1CP 

527 

8 oo 

-25.45 

8.98 x 1CP 

627 

900 

-28.90 

1.04 x id 7 

727 

1000 

-32.10 

l.o4 x id 7 

827 

1100 

-35.55 

1.16 x id 7 

927 

1200 

-39.39 

1.49 x id 7 

1027 

1300 

-42.30 

1.29 x id 7 

1127 

1400 

-45.90 

1.46 x 10 7 

1227 

1500 

-49.15 

1.45 x 10 7 

1327 

1600 

-53.25 

1.88 x 10 7 

1427 

1700 

-56.45 

1.81 x 10 7 

1527 

1800 

-60.60 

2.28 x 10 7 

1627 

1900 

-64.50 

2.63 x 10 7 

1727 

2000 

-67.40 

2.32 x id 7 

*TiCl 4 + 2 Ha 

+ 2 CO2 ^ TiOg 

+ 4 HCl + 2 CO 
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C. THERMODYNAMICS OF ZINC ORTHOTITANATE REACTIONS 

The suitability of the thermodynamics for the zinc-titania system 
has been established by thermodynamic calculation and previous investi- 
gations. Two reactions were considered for production of zinc ortho- 
titanate ; i . e . , 

2 Zn + TiCl 4 + 2 H 2 + 4 C0 2 + 4H, * Zn 2 Ti0 4 + 4 HC1,, + 4 TO (6) 
g g g g g s 8 S ' ' 

and 

2 ZnCl 2 + TiCl 4 + 4 C0 2 - Zn 2 Ti0 4 + 8HC1 + 4 CO . (7) 

egg S g g 

By summing the free energies of the products of the reactions and 
subtracting the sum of the free energies of the reactants, the free 
energy changes for Eqs. (6) and (7) were obtained from room temperature 
to 2000°K. These values are summarized in Tables III and IV. Free- 
energy data 9 5 10 were available for all the compounds in Eqs. (6) and 
(7) except Zn 2 Ti0 4 . An estimate for Zn 2 Ti0 4 was made using the follow- 
ing reaction: 


2 ZnO + Ti0 2 ^ Zn 2 Ti0 4 (8) 

This reaction occurs above 780°C 11 and thus the free energy of one 
mole of Zn 2 Ti0 4 must be more negative than the sum of the free energies 
of two moles of ZnO and one of Ti0 2 . Assuming that the free energy 
change for reaction (8) is at least -5,000 cal/mole Zn 2 Ti0 4 , a free- 
energy value for Zn 2 Ti0 4 was calculated, as summarized in Table V. 

As indicated in Tables III and IV, the free- energy changes for 
Eqs. (6) and (7) are substantially negative with large equilibrium 
constants above 1300°K (1023°C). Thus the thermodynamics are favorable 
for the formation of Zn 2 Ti0 4 above 1000°C where Zn 2 Ti0 4 is the only 
stable zinc titanate compound. 

Growth conditions for zinc orthotitanate within pressure and tem- 
perature ranges of the system were calculated by a CDC 6600 computer 
using program EQUICA. Program EQUICA is designed to compute the 
equilibrium compositions of a given set of reaction species by mini- 
mizing the Gibb’s free energy. The input data consisted of: 

a) possible species formed within the system 

b) Gibb's free energy values for these species 

c) temperature 
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Table III - Free -Energy Changes and Equilibrium Constants as Functions 
of Temperature for Zinc Orthotitanate Production Using 
Zinc Metal 


T 

(°c) 

T 

(°K) 

AG° 

Reaction* 

(kcal/mole) 

Reaction* 

27 

300 

- 38.90 

2.19 x 10 28 

127 

400 

- 46.25 

1.87 x 10 2S 

227 

500 

- 50.15 

8.36 x 10 21 

327 

600 

- 50.90 

3-48 x 10 18 

427 

700 

- 54.15 

8.09 x 10 16 

527 

800 

- 56.85 

3.40 x 10 15 

627 

900 

-59.50 

2.82 x 10 14 

727 

1000 

-61.80 

3.22 x 10 13 

827 

1100 

- 63.75 

4.65 x 10 12 

927 

1200 

- 66.19 

1.14 x 10 12 

1027 

1300 

-68.00 

2.71 x 10 11 

1127 

1400 

-69.70 

7.61 x IO 10 

1227 

1500 

-71.65 

2.76 x IO 10 

1327 

1600 

- 74.25 

1.39 x IO 10 

1427 

1700 

- 75.85 

5.65 x IO 9 

1527 

1800 

-78.30 

3.22 x IO 9 

1627 

1900 

-80.80 

1.97 x IO 9 

1727 

2000 

-81.90 

8.92 x IO 8 

* 

2 Zn + 

TiCl 4 + 2 H 2 + 4 C 0 2 

^TZn 2 Ti0 4 + 4 HC 1 + 4 CO 
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Table IV - Free-Energy Changes and Equilibrium Constants as Functions 
of Temperature for Zinc Orthotitanate Production Using 
Zinc Chloride 


T 

(°c) 

T 

(°K) 

AG° 

Reaction* 

(kcal/mole) 

K P 

Reaction* 

27 

200 

47.00 

5.27 x 10" 35 

127 

4oo 

31.55 

5.76 x 10“ 18 

227 

500 

16.45 

6.44 x 10” 8 

327 

600 

13.10 

I.69 X 10"5 

427 

700 

4.65 

3.53 x 10-2 

527 

800 

- 3.05 

6.81 x 10 

627 

900 

-10.70 

3.97 x 10 2 

727 

1000 

-17.80 

7.77 X 10 3 

827 

1100 

-26.91 

2.22 x 10 s 

927 

1200 

- 36.91 

5.28 x 1 CP 

1027 

1300 

- 45.28 

4.10 x 10 7 

1127 

1400 

- 54.54 

3.27 x 10 8 

1227 

1500 

- 63.25 

1.65 x 10 9 

1327 

1600 

- 70.85 

4.77 x 10 9 

1427 

1700 

- 76.85 

7.60 x 10 9 

1527 

1800 

- 84.30 

1.72 x IO 10 

1627 

1900 

-91.60 

3.45 x IO 10 

1727 

2000 

- 96.90 

3.89 x IO 10 

*2 ZnCl 2 

+ TICI4 + 4 h 2 + 4 

C 0 2 ^ Zn 2 Ti 0 4 + 8 HC 1 + 4 CO 
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Table V - Estimate of Maximum Free Energy of Zinc Orthotitanate 
(Assuming aGr = -5000 cal/mole Z 2 Ti04) 


Temperature 
°K °C 

2 AG° Zn0 
(kcal/mole) 

AG°Ti0 2 
( kcal/mole ) 

AG °Zn 2 Ti0 4 

Estimate* 

(kcal/mole 

300 

27 

-152.20 

-212.00 

-369.60 

400 

127 

-149.30 

-207.90 

-362.20 

500 

227 

-142.60 

-203.60 

-351.20 

6oo 

327 

-137.90 

-199.30 

-342.20 

700 

427 

-133.20 

-194.30 

-333.15 

8oo 

527 

-121.00 

-190.75 

-323.75 

900 

627 

-123.00 

-186.55 

314.50 

1000 

727 

-118.00 

-182.35 

-305.35 

1100 

827 

-112.20 

-178.00 

-295.20 

1200 

927 

-106.80 

-173.89 

-285.69 

1300 

1027 

-101.40 

-169.55 

-275.95 

l4oo 

1127 

- 95.80 

-165.45 

-266.25 

1500 

1227 

' - 90.20 

-161.20 

-256.40 

1600 

1327 

- 84.40 

-157.25 

-2*46.65 

1700 

1427 

- 78.80 

-153.10 

-236.90 

1800 

1527 

- 73.20 

-149.35 

-227.55 

1900 

1627 

- 67.60 

-145.25 

-217.85 

2000 

1727 

- 62.00 

-141.15 

-208.15 


AG °Zn 2 Ti0 4 “ 2 AG ZnO + AG Ti0 2 " 5000 cal / mole Zn 2 Ti0 4 
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d) pressure 

e) initial estimate of the composition of these species. 

The output consisted of: 

a) the equilibrium concentrations of the species 

b) equilibrium constants 

c) partial pressure of the species 

d) the total pressure. 

Figure 2 shows the equilibrium concentration of zinc orthotitanate 
for various temperatures and system pressures . Figure 3 shows the 
equilibrium constant for zinc orthotitanate formation for various 
temperatures . 


D. THERMODYNAMICS OF CALCIUM TUNGSTATE (CaW0 4 ) REACTIONS 
The computer output for the summary reaction 

CaCl 2 + WC1 5 + 1 H 2 + 4C0 2 * CaW0 4 + 7HC1 + 4C0 (9) 

shows thermodynamic feasibility at all temperatures, Table VI. Calcium 
chloride cannot be obtained by chlorinating calcium metal because the 
metal has a higher vapor pressure than CaCl 2 . Chlorination of CaO 
offers two possibilities, 

CaO + Cl 2 ^ CaCl 2 + l/2 0 2 (10) 

CaO + CO + Cl 2 ^ CaCl 2 + C0 2 (ll) 

Both reactions are thermodynamically favored at all temperatures, 

Tables VII and VIII. The vapor pressure of calcium chloride. Figure 4, 
requires temperatures above l600°C for sufficient volatility. 

Generating tungsten chloride, WC1 5 , by chlorination of tungsten 
metal is feasible due to the low vapor pressure of the metal. Various 
lower chlorides may be formed as intermediate compounds, but WC1 5 has 
the lowest free energy above 400°C, Figure 5. 

Growth conditions for CaW0 4 were calculated using program EQUICA. 
The equilibrium concentration of CaW0 4 as a function of temperature and 
pressure is shown in Figure 6. Figure 7 shows the equilibrium constant 
for CaW0 4 formation at various temperatures. 
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V P = 0.055 otm 
□ P = 0.110 otm 
A P = 0.220 atm 
o P = 0.440otm 


a 


B 


A 



O 


a 

v 


EQUILIBRIUM CONSTANT 



TEMPERATURE (°K) 


Fig* 3 - Equilibrium Constant as a Function of Temperature for 
ZngTi0 4 Formation 
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Table VI - Thermodynamic Data for 

Calcium 












Tungstate System (CaWO-i) 





THE i 

FOLLOWING DATA IS FOR 

THE 

REACTION 

■— 

- 





— 



reactants 

2 C ACL 2 

2 

WC L 5 

7 

H2 


8 CO? 





PRODUCTS 

2 CAW04 EST14 

HCL 

8 

CO 






. 


DELTA 

G 

IS 

-11200 

CALORIES 

AT 

300 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0. 14A5F + 09 


DELTA 

G 

IS 

-36000 

CALOR I5S 

AT 

400 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0 • 46 8 9 F -fRO 


DELTA 

G 

I S 

-70400 

CALORIES 

AT 

500 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0. 5946E *-31 


DELTA 

G 

IS 

-80200 

CALORIES 

AT 

600 

DEGREES 

KELVIN 

AND 

the 

EQUILIBRIUM 

CONSTANT 

IS 

0. 1641 E + 30 


DELTA 

G 

IS 

-95300 

CALOR IES 

AT 

700 

DEGREES 

KFLV IN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0.5708F+30 

H 

OELTA 

G 

is 

-110000 

calories 

AT 

800 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0. 11 30E+31 

-o 

DELTA 

G 

IS 

-1 24109 

CALOR IES 

AT 

900 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0. I374E.01 


OELTA 

G 

IS 

-1 37600 

CALORIES 

AT 

1000 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONST ANT 

IS 

0. 1 1 RSE*31 


DELTA 

G 

IS 

-152100 

CALORIES 

AT 

1100 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0. I667F + 31 


DELTA 

G 

IS 

-165500 

CALOR IES 

AT 

1200 

DEGREES 

KELVIN 

AND 

THE 

EQUI L I BRI UM 

CONSTANT 

IS 

0. 1394F. + 31 


DELTA 

G 

IS 

-176800 

CALORIES 

AT 

1300 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

! S 

0. 531 2 c ♦ 30 


DELTA 

G 

IS 

-18° 100 

CALORIES 

AT 

1400 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONST ANT 

IS 

0.3329F«-30 


DELTA 

G 

IS 

-200900 

CALOR IES 

AT 

1500 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0. 1R77P4-30 


DELTA 

G 

I s 

-216400 

CALORIES 

AT 

1600 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0. 3642E* 30 


OELTA 

G 

IS 

-226300 

CALOR IES 

AT 

1700 

DEGREES 

KELVIN 

ANO 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0.22 5 IE* 30 


. DELTA 

G 

IS 

-293000 

CALORIES 

AT 

1600 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0.3211E+30 


OELTA 

G 

IS 

-258900 

CALORIES 

AT 

1900 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0. 6064EO0 


OELTA 

G 

IS 

-271500 

CALOR IES 

AT 

2000 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0.4684E+30 



Table VII - Thermodynamic Data for Calcium 
Oxide Chlorination 


THE FOLLOWING DATA IS I OR THE REACTION 


R F. AC T AN 

TS 


pcnpiJCTS 


DELTA G 

IS 

- 

DELTA G 

I S 

- 

0 F I T *. G 

I S 

- 

0 FI T \ G 

IS 


0 F L T A G 

IS 

— 

DELTA r. 

I S 

- 

delta g 

IS 

- 

DFLTA G 

IS 

- 

OCLTA C, 

I S 

— 

DELTA G 

I S 

- 

0 r L T A G 

IS 


0 ( L T A G 

I S 

— 

DELTA G 

IS 

- 

DELTA G 

I s 

— 

DIET A G 

IS 

- 

'delta g 

IS 

- 

OFt TA G 

IS 

— 

DELTA G 

IS 

- 


1 C AQ 


1 Cl 2 


1 CACL2 


• 321 00 CAL OR I ES 
JJP S 0 _r A l CM' I r. s 
TO 3. TO CALORIES 
70 *-,00 CAl.HklES 

??'sob‘ ulni! ic‘s 

2 8- 1 0 9 f Al.nr IFS 
27900 CALORICS 
278 00 JCAI.OP I c_S 
27*59 CAL Of* IKS 
?7ono CALORIES 
? 7 SO 0 CA L OR i £ S 
27500 CALORIES 
27 000 CAL OP IE S 
27200 CALORIES 


-27 000 CALORIES’ 


I 02 

AT 300 

AT A 0 0 

A T 500 

AT 600 

_AT 7 00 

AT ROD* 
A T_ 000 
AT T 6 CO- 
AT 1100 
AT 1200 ' 

_ AT _noo_ 

AT VACO" 
AT_ I 500_ 
AT 1 6 00 
_AT 1700 
AT r?. 00 

_A T 19 00- 
AT 2000 


DEGREES 
DEGREES 
0 FORTES 
OGGi'L C S 
DEGREE S 

0 force’s 

DEGREES 
0 EGO E FS 
DEGREES 
OFGRC ES 

DEGREES. 

"nEGKTES 

DEGREES 

l) EG 2 FES' 

DEGREES 

DEGREES 

DEGRFJES 

DEGREES 


'<[- LV ID A ND THE ECU I LI 8 R I UM 

KELVIN A YD THE E QU I L I OR I U M 
KELVIN AND THE COD I L I RR I »JM 
KELVIN AMD T HE EQUILIBRIUM 

_K EJ-_V 1 : 1 AND T HA G OU 1 1. I B P I U 9 

KELVIN AND THF n l.'O I L I BR I UM 
KELV IJN_ANO THE EQUILIBRIUM 
KELvi N AMO THE " E QU I L'l RR I UM 
. K CIVIN_AN0 TH C - EQUILIBRIUM 
KELVIN AND THE EQUILIBRIUM 

KELV IN AND THE EQUILIBRIUM 

KFLVIN AND THE EOUI LIBRIUM* 

KELVIN _ A NO THF ECUI LI OR I IJ.M 

KFL VI N "AND "THE "cQUIL I BR I UM 

KFLVIN A ND-THF FOU I L I BR I UM 

KELVIN AND THE ECU I LI BR I UM 
KFLVIN A ND THE S OU I L I_B R I U M 
KELVIN AND THE ECU I LI RrTum 


CONSTANT IS 
CONSTANT Ts" 
CONSTANT IS 
CONSTANT" I S 

_C O N STAN T I S- 

CONST AN T IS 
CONSTANT^ IS 
CONSTANT I S ' 
CONSTANT IS 

CONSTANT* IS* 
C 0 N S T A N T_ I S 
CONSTANT ”l S 
-CONSTANT IS_ 
CONSTANT * I S 

JT. 0 N S T A N T I S 

CON St ANT IS*" 
_C ON ST ANT- IS 
constant is 


_0. 5226F 
0.*9295E 
_0. 2996F 
0. 99 36E 
_0. 5721E 
0". l 399F 
_0. 1959F 
0. 1 6 9 5 F* 
0. 3832F 
0. 12078* 
_0 • < 1 7 2 I E 
0. 2228E 
-0. 1 1 62E 
0. 5710E 
_0. 3933E 
0. 2065E 
.0. I 3 A 6 F 
0. 8926E 


26 ■ 
19i 
1 5 ■ 

12 i 

101 

09 1 
03 : 
07! 
06 ! 
06 i 
05 
05 ; 
05 i 

09"; 

09 I 
09 ! 
09 ! 

03" 



H 

VO 


Table VIII - Thermodynamic Data for Calcium 
Oxide Chlorination 



IMF FOLLOWING DATA IS FOR'thE REACT I CM 


R FACTAMTS 


1 CAC 


1 Cl 2 


1 CO 


PRODUCTS 1 C ACL? 

0 E L.TA . 0 _i s - 9 6 7 50 CALORICS 

DEI TA g IS - 93 A 00 'CAL>.P: ICS 

D F L M_G_. I S —90 409 C A I. OP I E S 

D F L TA 0 IS -8 7 2 0 6 "’ CA LORIES 

0 6 '. T_A_ G _I S - 84290 _CAt.Op [ cj 

OCLTA g" I s’’ -p'iibo CA topics' 

OF'.TA 1 } 1 S -75200 jCAl.op i r S 

0 Ft. T \ G I s -75 1 50 " C AL O R i F S 

OFt.TA _c, I s - 72 -S 0 C_CAt. OR ICS 

0 61. 1 A G IS -70 A OO CA LOR IF S 

o .T.L T _\ _0 .I S -6 81 5 0 CAL 0 R I F S_ 

DFl.TA G is -66350 CALORIES 

Fl.T A r, is < 1 ?30 C AI.ORI c S_ 

DELTA g IS -61700 CALORICS 

.0 F L T A G I S - 59 70 Q.CAin !K I C S 

OFLTA G IS -57550 CALORIES 
PFI TA._G. 1 S_ _-55350 CALORICS 
DELTA G IS -53850 CALORICS" 


1 CO? 

AT 300 

AT A 0 0 
A T _ 500 
AT" 6 00 
AT 7 00 
AT 5 00 

AT 900 

AT ion' 6 " 
AT l 1 00 
AT 1?00 

A I I 300 

AT 1 AGO” 
A T 1 5 00 _ 
AT 1600 
A T 1 7 00 _ 
AT 1800 
AT. 1°00 
A T' 2000' 


DEGREES 
DFG \ F F S 

0 FG-'FFS 
OF GREF S 
DEGREES 
"6 FOR F E s’ 
n FC.R PCs 

degree s' 

. 0 1 'C.R F C S 
DEG'-' EES 
n f G 8 F F S 
O FGRE E S 
OTG'i E F S_ 
DCGRFES 
OFC.RFF S_. 
D FOP E F S 
DFGRFES 
DEGREES’ 


JSELVI N 
K.C 1 .VI 0 
K c L V I N 
KELVIN 
KFLVI IV_ 

kelvif: 

K F IV I N 
K EL VI N 
KFL Vlft. 
KELVIN 
.KELVIN' 
KELVIN' 
KELV I 0 
KFLVI N 
KfLVt N. 
KELVIN' 
KELVIN 
KELVI F 


.A i 1 ID 
AND 
AND 
AMD 
AND 
AMD 
.AND. 
AND 
AND. 
AND 
AND 
•A NO 
AND 
AND 
A NQ_ 
AMD 

A. MP- 
A’ ND 


JHE 

THE 

THE 

THE 

THE 

THE 

THF 

THE 

THF 

THE 

THE_ 

THE 

J.HE. 

THF 

JI[F 

THE 

THE 

THE 


-JECU I I I DR HIM 
FOll I L IBP I UM 
...EOUIL I DR l UN 
ECU I L I DR I UM 
,-POUIL I 3RIUW 
EUUI L I OR I UM 
_F.pt! l.L I OK I (JM 
Entm I DR I G'M 
-ECU It I 0 R I U M 
ECU IL I BRIIJM 
-.F.CU I-L l.BR I UM 

FOUR I BR I uV~ 

_ F-CGUJ L_I. BRJ UM 
EOUI L I P.R I UM~ 
FP UIL1BRIUM_ 
EOUI LIORIUm’ 
..EOUIL I OR I UM 
E GU I L I BR I UM 


-COp.'STANT I S_ 

CONSTANT IS 
. CONST ANT_I S 
CONSTANT is" 
CONST AN T_ I S 
CONSTANT I S 
CONSJ ANT_I_S 
CONSTANT is" 
CONST AM T_.I S_ 
CONSTANT IS 

_C . P N S.T A n T I S_ 

CONSTANT IS 

CGNISJ ANT IS_ 

CONSTANT IS 
C.Q.NS TAN T_I S_ 
CONSTANT IS 
CONSTANT I S 
CONSTANT" I s’ 


_p. 30 . 7 _ 7 E 
0 . 10 S 5 E 
_ 0 . 3289 E 
0 . 5826 E 
- 0 . 1953 E 
0 . 1637 F 
— 0.97 9 8 F 
0 . 2663 E 
_ 0 . 2 66 A F 
0 . 6647 E 
0.287 IE 
0 . 2283 E 
- 0.-2 30 1 F 
0 . 2682 E 
_ 0_..47 3 8 E 
0.97 30 E 
. 0.2 3 29 F 
0 . 6597 F 


_7. l_ 
' 52 , 
_* 0 J 

■32{ 
-2 7 J 

■23, 
_1 9 ' 
i 7 ": 
..I 5 .; 

13 

-12 i 

HI 

_10J 
09 i 
__p 8 j 
07 

-P7.1 
Q 6 | 



Pressure 



Temperature, °K x 10" 2 


Fig. 4 - Vapor Pressure of CaCl 2 
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AG (kcal/mole) 



Fig. 5 - Free Energies of WCI25 WCI45 WC 1 5 , and WC 1 6 


EQUILIBRIUM CONCENTRATION (MOLES) 



TEMPERATURE (°K) 

Fig. 6 - Equilibrium Concentration of CaWO 

4 
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E. THERMODYNAMICS OF LANTHANUM OXIDE (La 2 0 3 ) REACTIONS 

The summary reaction for vapor phase production of lanthana powders, 

2LaCl 3 + 3H 2 + 3C0 2 * La 2 0 3 + 3C0 + 6HC1 (12) 

was introduced into the I EM 360/75 program for computing thermodynamic 
reaction free energies and equilibrium constants. The data output, 

Table IX, shows that temperatures above l600°C are necessary to make 
the reaction thermodynamically favored. Lanthanum chloride, LaCl 3 , 
can be obtained by chlorination of La 2 0 3 , 


La 2 0 3 + 3C1 2 ±; 2LaCl 3 + ^ 0 2 (13) 

which is thermodynamically feasible at all temperatures. Table X. The 
vapor pressure of lanthanum chloride, Figure 8, indicates that chlorin- 
ation above 1500°C will be required to ensure volatilization of the 
chloride. These restrictions on chlorination and reaction pose severe 
technical problems in the vapor phase reaction system now employed. 

Growth conditions for La 2 0 3 were calculated using program EQUICA. 
The equilibrium concentration of La 2 0 3 as a function of temperature and 
pressure is shown in Figure 9- Figure 10 shows the equilibrium constant 
for Ca 2 0 3 formation at various temperatures . 

In addition to the production of lanthanum chloride by reaction 
(l3)j lanthanum chloride was produced by the following reactions: 

La 2 0 3 + 3C1 2 + 3C ^2 LaCl 3 + 3 CO (l4) 

La 2 0 3 + 3C1 2 + 3C0 * 2LaCl 3 + 3C0 2 (15) 

Thermodynamic data for reactions (l4) and (15) are summarized in Tables 
XI and XII, respectively. 


LaOCl Reactions 


Since LaOCl occurred frequently in the La 2 0 3 runs, the thermodynamics 
of possible reactions were considered. Free-energy data were not avail- 
able for LaOCl and an estimate was made using the following reaction 

La 2 0 3 + LaCl 3 * 3LaOCl (l6) 

Assuming that the free-energy change is at least -10,000 cal/mole, a 
free-energy value for LaOCl was calculated and summarized in Table XIII. 
Possible reactions occurring in the chlorinator are: 
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Table IX - Thermodynamic Data for Lanthanum 
Oxide System (La 2 0 3 ) 


THE FOLLOWING DATA IS FOX THE REACTION 


REACTANTS 


2 L A CL 3 


3 . C02 


3 H2. 


PRODUCTS 
DELTA G IS 
DELTA 
DELTA 
DELTA 
OELTA 

delta 

DELTA 
DELTA 
.DELTA 
DELTA 
DELTA 
DELTA 
DELTA 
DELTA 
DELTA 
OELTA 
DELTA 
DELTA 


G 

G 

G 

G 

G 

G 

G 

G 

‘g 

G 

G 

G 

G 


IS 

I s 
i s 
is 
is 
is 

IS 

I s_ 
is 
is 
is 
i s_ 

IS 
I S 
IS 
IS 
IS 


1 LA 20 3 

121950 CALORIES 
"113100 CALORIES 

9S8pO_C ALOR I FS 

93700 CALORIES 
_ 8 3 7 50 CALORIES 
"72500 CALORIES 
_ 65 A 00 CALORIES 
55350 CALORIES 
*3 5 300 JCA 1081 £ S , 
3 7 A 00 CALORIES 
_ 316 50 CALORIES 
" ’ ? 6 1 0 0 CALORIES 

20150 CALORIES 

13100 CALORIES 

3000 CALORIES 

1650 CALORIES 
_ -6050 CALORIES 
-10250 CALORIES 


6 HCL 
AT 3 00 
AT" 600 

AT 5 00 

AT 6 00 
AT 700 


AT 

AT 

AT 

A L 

AT 

AT 

AT 


800 
900 
1000 
.11 00 
1200 
13 00 
1600 


AT 1500 
AT 1600 
AT 1700 
AT" 1800 
A T _1 900 
AT 2000 


OEGREES 

DEGREES 

DEGREES 

DEGFEES 

DEGREES 

DEGREES 

degrees 

DEGREES 

DEGREES 

'DEGREES 

OEGREES 

'DEGREES 

DEGREES 

DEGREES 

DEGREES 

'degrees 

OEGREES 

OEGREES 


CO 

KELVIN 
KfcLVI N 
KELVIN . 
KEL VI N 
KELVIN, 
KELVIN 
KELVIN- 
KELVIN 
KELVIN. 
‘ KELVIN 
KF L VI N 
"KELVIN 
KELVIN 
KFLVI N 
KELVIN 
KELVIN 
KELVIN 
KELVIN 


AND 
AND 
AND 
A NO 
AND 
AND 
AND 
AND 
AND 
AND 
AND 
AND 
AND 
AND 
.AND 
AND 
AND 
"A NO 


THE 

the 

THE_ 

THE 

THE. 

THE 

THE. 

THE 

THE 

"THE 

THE 

"the" 

THE 

THE 

JHE 

the 
_t h e_ 
THE 


CCUILI 
EOUILI 
FGUIl I 
EUUILI 
EOUILI 
ECU I L I 
EOU 1 1. 1 
’ ECUIL 1 
EOU I L I 
" EOUIL I 
EOUILI 
EOUILI 
EOUILI 
FQUIL I 
JEOLIILJ. 

EOUILI 
_£CU l L 1 
EOUILI 


(3R IUM 
BRI UM 
13RI UM. 
BRI U M 
BR 1 UM. 
BRI UM 
P R I UM 
BF l UM 
BRI UM 
BRI UM 
DR IUM 
BRIUM 
BR I UM 
BP I UM 
BRI UM 
BRIUM 
BRIUM 
BRIUM 


CONST ANT_I S . 
CONSTANT is 
CONSTANT... IS.. 
CONSTANT IS 
CONSTANT _ I S 
CONSTANT IS 
CONSTANT .IS 
CONSTANT IS 
CONSTANT. I S_ 
CONSTANT IS 
CONST ANT I S. 
CONSTANT IS 
CONSTANT...! S. 
CONSTANT IS 

CON STAN T I.S_ 

CONSTANT IS 
CONSTANT IS_ 
CONSTANT IS 


0.0 

0. 158 r -E-61 
0. 66 716^3. 
0. 73 62 E -36 
0. 70 78F-2 6 . 
0. 1557E-19 
0. 1310E-15. 
0.7985E-12 
_0 . 99 77 F — 09. 
0. 1562E-06 
0.6772E-05. 
0. 861 9E-06 
_0. 11 59E-02. 

0. 1626E-01 
_0. .93 66f -01. 
0.6306E 00 
0.36136 01 
0.1319E 02 



Table X - Thermodynamic Data for Lanthanum 
Oxide Chlorination 


THE FOLLOWING DATA IS FOR THE REACTION 
•‘RE" ACTANTS 1 ."LA 203 -3 CL-2 

"PRODUCTS 2 L ACL 3 ~ i,02 


DELTA 

G 

IS 

-'79000 

CALORIES 

A* 

300 

DEGREES 

KELVIN 

AND 

the 

EQUILIBRIUM 

CONSTANT 

IS 

0.8193E 

59 

DELTA 

G 

IS 

- 7 2 0 0 0 

CALORIES 

AT 

900 

DEGREES 

KELVIN 

AND 

THE 

FCUIL l OPIUM 

CONSTANT 

IS 

0.6G29E 

90 

DELTA 

G 

IS 

-70300 

CALORIES 

AT 

500 

DEGREES 

KELVI N 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0.5378E 

31 

DELTA 

G 

IS 

-66800 

CALORIES 

AT 

600 

DEGREES 

KFLVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0. 11 55F 

26 

DELTA 

G 

I s 

— 60 500 

CALORIES 

AT 

700 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0.S807F 

21 

DELTA 

’ G 

IS 

-6260 0 

‘CALORIES 

AT 

800 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0. 1267F 

IB 

DELTA 

G 

IS 

-62700 

CALORIES 

AT 

900 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CON SI ANT 

IS 

0- 1 6 8 6 F 

16 

DELTA 

G 

IS 

-59900 

CALORIES 

AT 

1000 

DEGREES 

KFLVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0.9615E 

13 

DELTA 

G 

IS 

-56700 

calories 

AT 

1100 

DEGREES 

KFLVIN 

AND 

THE 

EOUIL IBP. IUM 

CONSTANT 

IS 

0. 1S96F 

12 

delta 

G 

IS 

-56300 

CALORIES 

AT 

1200 

DEGREES 

KFLVIN 

AND 

THE 

EQUILIBRIUM 

CONST ANT 

IS 

0. 1797E 

1 1 

DELTA 

G 

I S 

-57000 

CALORIES 

AT 

1300 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

I s 

0. 3831F 

10 

DELTA 

G 

IS 

-53500 

CALORIES 

AT 

1900 

OcGPEFS 

KELVIN 

AND 

THE 

ECU I L IBRIUM 

CONSTANT 

I s 

0.13S8E 

io 

DELTA 

G 

IS 

-59300 

CALORIES 

AT 

1500 

OEG°cES 

KELVI N 

AND 

THE 

FCUIL I8R I UM 

CONSTANT 

I s 

0.9372E 

09 

DELTA 

G 

IS 

—59900 

CALORIES 

AT 

1600 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0- 1523E 

09 

DELTA 

G 

I S 

— 61 A 00 

CAL OR I E S 

AT 

1700 

DEGREES 

KELVI N 

AND 

THE 

EQUILIBRIUM 

constant 

IS 

0. 7837E 

08 

DELTA 

G 

IS 

-62900 

CALORIES 

AT 

1800 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0. 3776E 

08 

OELTA 

G 

IS 

-63900 

CALORIES 

AT 

1900 

DEGPETtS 

KFLVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

_0. 1«R6 5F 

08 

DELTA 

G 

IS 

-69000 

CALORIES 

AT 

2000 

DEGREES 

KELVIN 

AND 

THE 

EQUILIBRIUM 

CONSTANT 

IS 

0.9867E 

07 



Pressure 



EQUILIBRIUM CONCENTRATION (MOLES) 



Fig. 9 - Equilibrium Concentration of La 2 0 3 as a Function of 
System Conditions 
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EQUILIBRIUM CONSTANT 



Fig. 10 - Equilibrium Constant as a Function of Temperature 
for La 2 0 3 
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Table XI - Free-Energy Change and Equilibrium Constant 
for Chlorination of La 2 C>3 with Carbon* 


T(°C) 

T(°K) 

ag° 

Reaction 

(kcal/mole) 

K P 

Reaction 

227 

500 

-I8I.6O 

7.24 x 10 75 

327 

600 

-186.85 

1.16 x 10 68 

427 

700 

-191.15 

4.84 x 10 59 

527 

800 

-193.70 

8.33 x IO 52 

627 

900 

-200.25 

4.28 x 10 48 

727 

1000 

-203.25 

2.65 x 10 44 

827 

1100 

-207.00 

1.35 x 10 43 

927 

1200 

-212.75 

5.63 x 10 38 

1027 

1300 

-200.05 

9-73 x 10 36 

1127 

i4oo 

-227.25 

3.01 x 10 35 

1227 

1500 

-234.50 

1.48 x 10 34 

1327 

1600 

-241.70 

1.04 x 10 33 

1427 

1700 

-249.35 

1.15 x 10 32 

1527 

1800 

-256.35 

1.34 x 10 31 

1627 

1900 

-263.65 

2.13 x 10 3 ° 

1727 

3000 

-270.25 

3-42 x 10 29 


* La 2 0 3 + 3C1 2 + 3C ±5 2LaCl 3 + 3C0 
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Table XII - Free -Energy Change and Equilibrium Constant 
for Chlorination of La203 with CO* 


T(°C) 

T(°K) 

AG° ' 
Reaction 
(kcal/mole) 

*P .. 
Reaction 

327 

600 

- 234.10 

7.24 x 10 75 

427 

700 

-225.35 

2.31 x 10 7 ° 

527 

800 ' 

-215.O6 

5.50 x 10 58 

627 

900 

-208.80 

5,10 x 10 5 ° 

727 

1000 

-199.35 

3.73 x 10 43 

827 

1100 

-190.20 

6.20 x 10 37 

927 

1200 

-183.80 

3-00 xlO 33 

1027 

1300 

-178.05 ' 

8.62 x 10 29 

1127 

1400 

- 174.06 

1.46 x 10 27 

1227 

1500 

- 168.35 r. 

3.39 x 10 24 

1327 

1600 . 

-162.50 

1.58 x 10 22 

1427 

1700 

-158.00 

2.06 x IO 20 

1527 

1800 

- 153.15 

3.95 x 10 18 

1627 

1900 

- 147.85 

1.02 x 10 17 

1727 

2000 

- 142.75 

3.98 x 10 15 

* La 2 0 3 + 3 C 1 2 + 3 CO 

* 2 LaCl 3 + 3 C 0 2 
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Table XIII - Estimate of Maximum Free-Energy of LaOCl* 


T(°K) 

T(°C) 

AG° LaCl 3 
(kcal/mole) 

AG° La 2 0 3 
(kcal/mole) 

AG° LaOCl 
Estimate 
(kcal/mole) 

300 

27 

-238.30 

-402.60 

-223.63 

400 

127 

-233.00 

-393.10 

-218.70 

500 

227 

-227.40 

-384.50 

-213.97 

600 

327 

-222.70 

-376.60 

-209.77 

700 

427 

-217.30 

-368.10 

-205.13 

800 

527 

-211.30 

-360.00 

-200.43 

900 

627 

-207.00 

-351.30 

-196.10 

1000 

727 

-200.90 

-342.40 

-191.10 

1100 

827 

-195.80 

-334.90 

-186.90 

1200 

927 

-191.60 

-326.90 

-182.83 

1300 

1027 

-188.00 

-319.00 

-179-00 

l4oo 

1127 

-184.50 

-310.50 

-175.00 

1500 

1227 

-180.90 

-302.50 

-171.13 

1600 

1327 

-177-20 

-294.50 

-167.23 

1700 

1427 

-173.70 

-286.00 

-163.23 

1800 

1527 

-170.20 

-278.00 

-159-40 

1900 

1627 

-166.70 

-270.00 

-155.57 

2000 

1727 

-163.00 

-262.00 

-151.67 

* AG° LaOCl 

= 1/3 AG° 

La 2 0 3 + l/3 AG° 

LaCl 3 = -10,000 ' 

cal/mole LaOCl 
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La 2 0 3 + Cl 2 + C ft 2LaOCl + CO (17) 

La 2 0 3 + Cl 2 + CO ft 2LaOCl + C0 2 (18) 

La 2 0 3 + Cl 2 ft 2 LaOCl + 1/2 0 2 . (19) 

Thermodynamic data for reactions (17), (l8) and (19) are summarized in 
Tables XIV, XV, and XVI, respectively. Possible reactions occurring in 
the reaction zone are: 

2LaOCl + H 2 + C0 2 ft La 2 0 3 + 2HC1 + CO (20) 

LaOCl + LaCl 3 + 2H 2 + 2C0 2 ft 2La 2 0 3 + 4HC1 + 2C0 (21 ) 

LaCl 3 + C0 2 + H 2 ft LaOCl + CO + 2HC1 . (22) 

Thermodynamic data are summarized in Tables XVII, XVIII and XIX. 


F. REACTOR DESIGN, MIXING AND HEAT TRANSFER 

Prior experience in the area of particle production by vapor-phase 
reaction has shown that mixing of reactants was limited. To improve 
the mixing of reactant gases, a detailed study was performed and the 
results were used to design a reactor to provide suitable mixing of the 
gaseous species . 


i. Thermodynamic and Transport Properties 

Energy balances, fluid-mechanical, mixing and heat transfer 
calculations require the knowledge of thermodynamic and transport 
properties . Densities of reactants and products as well as heat capac- 
ities, viscosities, thermal conductivities and mass diffusivities were 
calculated or estimated whenever not available in literature. 

a. Viscosity 

For the estimation of gas viscosity at low densities, the 
following formula was used , 12 

4 = 2.6693 x 10” 5 ^5L ( 23 ) 

^4 
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Table XIV - Free-Energy Changes and Equilibrium Constants 
for LaOCl* 


T(°C) 

T(°K) 

£0° 

Reaction 

(kcal/mole) 

*P 

Reaction 

27 

300 

-77.47 

2.75 x 10 56 

127 

400 

-79-30 

2.14 x 10 48 

227 

500 

-80.53 

1.60 x 10 35 

327 

600 

-82.28 

9.43 x 10 29 

427 

700 

-83.72 

1.38 x 10 26 

527 

800 

-84.57 

1.27 x 10 23 

627 

900 

-86.75 

1.17 x 10 21 

727 

1000 

-87.75 

1.51 x 10 19 

827 

1100 

-89.OO 

4.83 x 10 17 

927 

1200 

-90.92 

3.63 x 10 16 

1027 

1300 

-93.35 

4.95 x 10 ls 

1127 

1400 

-95.75 

8.88 x 12 14 

1227 

1500 

-98.17 

2.01 x 10 14 

1327 

1600 

-100.57 

5.47 x 10 13 

1427 

1700 

-103.12 

1.8l x 10 13 

1527 

1800 

-105.45 

6.37 x 10 12 

1627 

1900 

-107.88 

2.57 x 10 12 

1727 

2000 

-no. 08 

1.07 x 10 12 

Xj 8*203 + 

CI2 + C i: 

2La0Cl + CO 
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Table XV - Free -Energy Changes and Equilibrium Constants 
for LaOCl* 


T(°C) 

t(°k) 

ag° 

Reaction 

(kcal/mole) 

Reaction 

27 

300 

-106.12 

7.24 x 10 75 

127 

400 

-103.60 

4.07 x 10 5S 

227 

500 

- 100.73 

1.08 x 10 44 

327 

600 

- 98.03 

5.15 x 10 35 

427 

700 

- 95.12 

5.00 x 10 28 

527 

800 

- 91.67 

1.11 x IO 25 

627 

900 

-89.6O 

5.75 x 10 21 

727 

1000 

- 86.45 

7.86 x 10 18 

827 

1100 

- 83.40 

3-73 x 10 16 

927 

1200 

- 81.27 

6.34 x 10 14 

1027 

1300 

- 79.35 

2.19 x 10 13 

1127 

1400 

-78.00 

1.51 x 10 12 

1227 

1500 

-76.12 

1.23 x 10 11 

1327 

1600 

- 74.17 

1.35 x IO 10 

1427 

1700 

- 72.67 

2.20 x 10 9 

1527 

1800 

- 71.05 

4.24 x 10 s 

1627 

1900 

-69.28 

9.34 x 10 7 

1727 

2000 

- 67.58 

2.43 x 10 7 

* L 3 

+ Cl 2 + CO 

- 2 LaOCl + C 0 2 
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Table XVI - Free-Energy Changes and. Equilibrium Constants 
for LaOCl* 


T(°C) 

t(°k) 

ag° 

Reaction 

(kcal/mole) 

K P 

Reaction 

27 

300 

- 44.67 

3.48 x 10 31 

127 

400 

- 44.30 

1.61 x 10 24 

227 

500 

- 43.43 

9.70 x 10 18 

327 

600 

- 42.93 

4.37 x 10 15 

427 

700 

- 42.16 

1.47 x 10 13 

527 

800 

- 40.87 

1.47 x 10 11 

627 

900 

- 40.90 

8.56 x 10 9 

727 

1000 

- 39*80 

5.00 x 10 s 

827 

1100 

-38.90 

5.36 x 10 7 

927 

1200 

- 38.77 

1.15 x 10 7 

1027 

1300 

- 39*00 

3.61 x 10 s 

1127 

1400 

- 39*50 

1.47 x 10 6 

1227 

1500 

- 39-77 

6.23 x 10 5 

1327 

1600 

- 39.97 

2.88 x 10 5 

1427 

1700 

- 40.47 

1.60 x 10 5 

1527 

1800 

- 40.80 

9.00 x 10 4 

1627 

1900 

- 41.13 

5-39 x 10 4 

1727 

2000 

- 41.33 

3.29 x 10 4 

^ L&2O3 + 

Cl 2 * 2 LaOCl 

+ 1/2 o 2 
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Table XVII - Free-Energy Changes and Equilibrium Constants 
for La203 Production from LaOCl* 


T(°C) 

t(°k) 

AG° 

Reaction 

(kcal/mole) 

K P 

Reaction 

27 

300 

60.62 

6.88 x IO -45 

127 

400 

57.70 

2.96 x 10" 32 

227 

500 

52.93 

7.26 x 10- 24 

327 

600 

51.23 

2.17 x 10“ 17 

427 

700 

47.92 

1.09 x 10 -15 

527 

800 

44.17 

8.58 x 10 -13 

627 

900 

41.80 

7.06 x 10- 11 

727 

1000 

38.45 

3.95 x 10 - 9 

827 

1100 

35.10 

1.06 x 10- 7 

927 

1200 

32.47 

1.22 x 10- 6 

1027 

1300 

30.55 

7.31 x 10- 6 

1127 

l4oo 

28.76 

3.31 x 10- 5 

1227 

1500 

26.72 

1.28 x 10 -4 

1327 

1600 

24.37 

4.69 x 10- 4 

1427 

1700 

22.67 

1.22 x 10- 3 

1527 

1800 

20.55 

3.20 x 10- 3 

1627 

1900 

18.38 

7.68 x 10- 3 

1727 

2000 

16.58 

1.54 x 10- 2 

* 2 LaOCl 

+ H2 + CO2 

i; La 203 + 2 HC 1 + 

CO 
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Table XVIII - Free-Energy Changes and Equilibrium Constants 
for La 2 03 Production* 


T(°C) 

T(°K) 

ag° 

Reaction 

(kcal/mole) 

*P 

Reaction 

27 

300 

91.23 

3.40 x 10- 67 

127 

400 

85.40 

2.17 x 10“ 47 

227 

500 

75.89 

6.85 x 10" 34 

327 

600 

72.47 

4.00 x 10 -27 

427 

700 

65.83 

2.78 x 10" 21 

527 

800 

58.33 

1.16 x 10- ls 

627 

900 

53.60 

9-62 x 10- 14 

727 

1000 

46.90 

5.61 x 10- 11 

827 

1100 

40.20 

1.03 x 10“ 8 

927 

1200 

34.93 

4.34 x 10- 7 

1027 

1300 

31.10 

5.90 x 10- 6 

1127 

1400 

27.40 

5.28 x 10" 5 

1227 

1500 

23.43 

3-85 x 10- 4 

1327 

1600 

18.73 

2.76 x 10" 3 

1427 

1700 

15.33 

1.07 x 10 2 

1527 

1800 

11.10 

4.49 x 10" 2 

1627 

1900 

6.77 

1.67 x 10- 1 

1727 

2000 

3.17 

4.51 x 10" 1 

* LaOCl + 

LaCl 3 + 2H-2 + 

2C0 2 La 2 0 3 + 4HC1 + 2C 
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Table XIX - Free-Energy Changes and Equilibrium Constants 
for LaOCl Production from LaCl 3 * 


T(°C) 

T(°K) 

ag° 

Reaction 

(kcal/mole) 

*P 

Reaction 

27 

300 

30.62 

4.94 x 10- 23 

127 

400 

27.70 

7.31 x 10" 16 

227 

500 

22.93 

9.44 x 10" 11 

327 

600 

21.23 

1.84 x 10" 8 

427 

700 

17.92 

2.54 x 10" s 

527 

800 

. 14.17 

1.35 x 10- 4 

627 

900 

11 .80 

1.36 x 10 -3 

727 

1000 

8.45 

1.42 x 10~ 2 

827 

1100 

5.10 

9.70 x 10- 2 

927 

1200 

2.47 

3.55 x 10' 1 

1027 

1300 

0.55 

8.08 x 10" 1 

1127 

i4oo 

-1.35 

1.60 

1227 

1500 

-3-28 

3-01 

1327 

1600 

-5.63 

5.88 

1427 

1700 

-7.33 

8.77 

1527 

1800 

-9.45 

1.40 x 10 1 

1627 

1900 

-11.62 

2.17 x 10 1 

1727 

2000 

-13.42 

2.93 x 10 1 

* LaCl 3 + 

C0 2 + h 2 

LaOCl + CO + 2HC1 
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in which n [= ]g cm -1 sec -1 , T[=]°K, ct[=]A and Hp is a slowly varying 
function of the dimensionless temperature KT/e . 13 Values of e/K and 
collision diameter, <x, were calculated from 


e/K = 1.15 Tb, a = 1.166 (V h )^ 3 

u liq 


( 24 ) 


in which e/K and T^ are in °K, cr[=]A., and V[=]cm 3 (g-mole) -1 . For 
multicomponent gas mixtures, the semi -empirical formula of Wilke 14 was 
used. 


^ix 


£ x i^i 

1=1 £ *0%d 

0=1 


(25) 


where 



and n is the number of chemical species in the mixture; x^ and xj are 
the mole fractions of species i and j; pj_ pj and M^, Mj are the vis- 
cosities and molecular weights respectively. 

b . Thermal Conductivity 

The Eucken approximation for estimating thermal conduc- 
tivity of a polyatomic gas at low density was used . 15 

k = (C p + 5/4 R) p/ M (27) 

where 

p = viscosity 
M = molecular weight 
R = gas content 
Cp = heat capacity. 
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For a gas mixture (n-ccanponents), 


n 


^mix = 






( 28 ) 


£ ViJ 


0=1 


where Xj_, Xj and $j_j are as in equations (25) and (26). 


c . Mass Diffusivity 

The following formula for the estimation of mass diffu- 


sivity of gases at low density was used 


is 


D AB 


0.0018583 



P °AB % 


where 


(29) 


TAB = V 2 (°A +0 b)> (30 ) 

T[=]°K, Ma, Mg are molecular weights of A and B, p[=]atm, ctab[=]A-> 
is a dimensionless function of temperature, and Dab[= 1 cm 2 / sec. A 
summary of transport properties are listed in Table XX. 


2. Reactor Gas Mixing 

Past experience in the area of gas phase particle production 17 
had shown that mixing of reactants was indeed limited. A detailed study 
of this problem was undertaken to design a reactor which would provide 
adequate mixing of the reactants and increase the product yield. The 
example product was CaW0 4 , produced by the overall reaction. 

CaCl 2 + WC1 5 + 7/2 H 2 + 4C0 2 * 4CO + CaW0 4 + 7HC1 (31) 

Thermodynamic calculations, previous experience and existing experimental 
facilities dictated the choice of operating conditions, as well as the 
reactor configuration. Fluid mechanical and mixing calculations dete rmin ed 
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Table XX - Estimated Transport Properties 


Species 

Viscosity 

h 

g cm" 1 sec -1 

Mass Diff. 
d AB 

cm 2 /sec 

Temperature 

°K 

Kinematic 
Viscosity 
p/p , cm 2 /sec 

CaCl 2 

3.49 x 1CT 4 


2,073 

5-33 x 10" 2 

WC1 5 

6.20 x 10" 4 


2,073 

2.91 x 10" 2 

CaCl 2 | 
WC1 5 ) 


(at p = 1.0 atm) 
0.37 

2,073 


Op 

8.0 x 10" 4 


2,073 


C0 2 

4.9 x 10“ 4 


2,073 

3.26 

co 2 

5.24 x 10 -4 


1,598 


H 2 

2.65 x 10" 4 


1,598 


ZnCl 2 

5.10 x 10" 4 


1,598 


TiCl 4 

4.60 x 10" 4 


1,598 



the efficiency for the degree of mixing under the chosen operating con- 
ditions and reactor configuration. 

In Figs. 11 and 12, the proposed reactor configuration is shown. 
One gas stream is composed of 0 2 + C0 2 , and the other gas stream is com- 
posed of CaCl 2 + WC1 5 + H 2 . Two heat sources were considered: an 

induction furnace and a oxy-hydrogen flame. A preliminary reactor design 
was based on stoichiometric reactant flow rates to give a first estimate 
of reactor dimensions. Heat balance calculations gave information of 
the additional flow rates of oxy-hydrogen flame. The oxy-hydrogen flame 
would provide heat for the reaction and create instability and turbulance 
at the central jet. Revised calculations based on the new rates provided 
the final reactor configurations and dimensions . 

Preliminary Reactor Design : The selected operating conditions 

were : 


Pressure: 0.5 atm 

Temperature: 2,074°K 
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Gas flow rates : 


C0 2 

CaCl 2 

WC1 5 

H 2 


800 cc/min (STP) 
200 " " " 

200 " " " 

700 " " " 


Stream 1 forms a jet issuing from a l/2" I.D. tube (central jet) with 
a volumetric flow rate of V x = 1.2 x 10 4 cc/min and an average issuing 
velocity of ^ = 1.57 x 10 4 cm/sec under the selected operating condi- 
tions. Stream 2 is split into eight jets issuing from eight l/4" orifices 
(peripheral jets) with a total volumetric flow rate of V 2 = 1.66 x 10 4 cc/ 
min and an average issuing velocity of 1.11 x 10 2 cm/sec. 


The density and kinematic viscosity of Stream 2 was calculated 
as P 2 = 2.18 x 10 -4 g/cc and v 2 = 1.7 cm 2 /sec respectively; and that of 
Stream 1 as P x = 1.2 x 10 -4 gm/cc and v x = 3.26 cm 2 /sec, respectively. 
Hence the Reynolds numbers of the jets are based on their issuing velo- 
city and diameter. 


Central jet: 


Peripheral jets: 


( NEe)l . (157) (l/2 6 * . 49 ( 32 ) 

(NRe )z = 2 : = 4l (33) 


Since the transition Reynolds number for a free jet is approximately 
30, turbulent velocity fluctuations should be present outside of the 
jets potential core. Assuming that Stream 1 has fully occupied the 
cross sectional area of the divergent nozzle of the reactor (Fig. 11); 
the time (equivalently, the length of the mixing zone) for a prescribed 
degree of mixing may be determined. 

The terminology of mixing may cause confusion, since there is 
generally a free interchange of words (e.g., blending, mixing, disper- 
sion) and their meanings. Brodkey 18 defines "mixing to mean any blending 
into one mass and mixture as a complex of two or more ingredients which 
do not have a fixed proportion to one another and which, however 
thoroughly comingled, are conceived as retaining a separate existence". 
The ultimate in any mixing process is submicroscopic homogeneity with 
molecules uniformly distributed over the field. Molecular diffusion 
will always be the controlling factor and when it is very slow, turbu- 
lence can speed up the diffusion by increasing the transfer area. It 
is the dispersive action inherent in turbulence that enhances the 
molecular diffusion. 


Criteria characterizing the quality of mixing are defined for 
the purpose of this report; 19 their mathematical treatments are avail- 
able elsewhere. 20-2 '’ 
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Scale of segregation : A measure of the size of the unmixed 

clumps of the pure components. As these clumps are dispersed, the 
scale of mixing is reduced. A linear scale, L s , is defined as 




g s (r) dr 


( 34 ) 


where , 


B s ( r ) 



(35) 


and 

a = c-c = concentration fluctuation 
c = instantaneous concentration 
c = time-average concentration 
a ' = rms fluctuation s 4W 2 . 


A volume scale, V s , is defined as 



where x~ and r~ are vector distances and r = ('r | . 

Intensity of segregation : The effect of molecular diffusion 
on the mixing process is a measure of the difference in concentration 
between neighboring clumps of fluid. The intensity of segregation is 
defined as 



(37) 


The intensity of segregation is unity for complete non -mixing and equal 
to zero when the mixture is perfectly uniform (a* 2 = 0). Figure l4 is 
a schematic representation of the scale and intensity of segregation 
during a gas mixing process. 
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The statistical theory of turbulence suggests the following 

relations 



( 38 ) 


where 


t = time of mixing 
t = time constant 


->4 (If 




In N 


sc 


] 


(39) 


and 

u = kinematic viscosity = p/p[= ]cm 2 /sec 
N sc = Schmidt number = -q/D 
D = molecular mass diffusivity [=]cm 2 /sec 
L s = scale of segregation 


The scale of segregation, L s , can be estimated. 



Af_ _ 0-3^1 r Q 
lOo " h’ 


A 3 = 15u u* 2 /e 


(40) 


(41) 


where 


A = a microscale of mixing 

u' = Vd 2 = rms velocity fluctuation in the axial direction 

In summary, the key to the mixing is the determination of the 
time constant, t. From equation (38), the time of mixing may be calcu- 
lated for any desired value of the intensity of segregation I g , i.e,, 
degree of mixing. 
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Mixing Calculations 


Values of the intensity of segregation were calculated for 
different values of the rms value of the axial velocity fluctuation 
u ' . The geometry is shown in Figs . 11-13 • 



Fig. 13 - Geometry Used for Mixing Calculations 



Fig. 14 - The Intensity and Scale of Segregation 
During Mixing 28 



Two cases were considered: 


(a) Mixing of the peripheral jets (l/4" nozzle) with 
surrounding gas considered as stagnant . Three values of u’ were 
assumed; u x = 0.0^V 2 , u 2 = 0.025v 2 and u 3 = 0.01v 2 , where v 2 = 110 cm/ 
sec. Corresponding time constants were calculated: t x = 0.236 sec, 
t 2 = 0.472 sec, and t 3 = 1.165. Figure 15 summarizes the results of 
this case where the intensity of segregation is plotted against time. 

(b) Mixing of the central jet with surrounding gas con - 

sidered stagnant . Time constants of mixing were calculated as in the 
preceeding case with v x = 157 cm/sec. The corresponding time constants 
were calculated as: t x = 0.110 sec, t 2 = 0.332 sec and t 3 = O.83O sec. 

Figure 16 summarizes the results of this case. 

The length of the reaction zone was calculated on the basis of 
mixing time (for a given intensity of segregation and an average flow 
velocity of both gas streams). For an 80f 0 degree of mixing, i.e., 

I s = 0.2, and an average gas velocity of 23.5 cm/sec the time of mixing 
was calculated as t = 1.875 sec. Taking N = 2*, the dimensionless 
measure of mixing distance, cr. 


cr = k 0 u't 


(42) 


drops from cr = 6 to cr 3 4.5*. Hence the actual mixing time is 


tmix - 1*875 f j = 1.40 sec 


(43) 


and the length of the mixing zone 


L,,! = (23.4 cm/sec) (l.4 sec) 3 32.8 cm 


(44) 


These preliminary reactor design data formed the basis for subsequent 
reactor design. 

Reactor Design 

The rms value of the axial velocity fluctuation, u', varies 
with distance downstream in the reactor. Moreover due to the spreading 
of the central and peripheral jets, their average velocity also varies 
with distance downstream in the reactor. The spreading of the jets is 
not free, due to interactions of two adjacent peripheral jets and the 


* [29] 
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Fig. 15 - Preliminary Reactor Design. Mixing of Stream 2 
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Fig. 16 - Preliminary Reactor Design. Mixing of Stream 1 
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central jet. Thus, the flow field becomes extremely complicated and 
almost hopeless from an exact analytical point of view. 

Semi-empirical approximations have been used as well as reason- 
able simplifications to make the problem tenable. An energy (heat) 
balance around the reactor established the heat requirements being 
supplied by the oxy-hydrogen flame, thus revising the flow rates in 
the central jet. 

There is a wealth of information concerning free turbulence 
and especially round free jets. 30-32 Figures 17 and 18 show the 
nomenclature and configurations used in these calculations . 

The jet issues with a velocity Vp from an orifice with diameter, 
d, into a stream with uniform velocity v s . In the jet region, the 
axial velocity component v x is assumed to be composed of the velocity 
v s of the ambient stream and the velocity V lf which is the deviation 
from the velocity v s caused by the jet (Fig. 17)- Close to the orifice, 
the turbulent mixing zone has zero width W; this width W(x) increases 
with increasing distance from the orifice to a distance x c where the 
mixing zone covers the entire jet (Fig. 18). If the velocity distribu- 
tion at the orifice is uniform, the velocity at the axis maintains its 
constant value vp up to this distance x c . Up to x c , the jet contains 
a central potential core region. Farther downstream the velocity at 
the axis of the jet decreases with increasing distance. 

The length of the core region, xc, is determined mainly by the 
value of the ratio p s v s /vpj one empirical correlation 33 is: 


-£=4 + 12m (45) 

d 

The velocity of the jet at the axis decreases hyperbolically with in- 
creasing downstream distance x, 


v i >max 



s 


x c 

— for x > X c 


(46) 


and the change in jet width, W, is given by 

| = (2)l-M for * > x c . 


(47) 
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The heat of reaction of the oxy-hydrogen flame was calculated 
as: -60.298 kcal/g-mole at T = 2,073°K. From the heat balance calcu- 

lations around the reactor it was determined that the amount of heat 
supplied by the oxy-hydrogen flame is 


Qf = 40 kcal/hr (48) 

Hence the volumetric flow rates were revised. 

V3. = 1.373 x 10 4 cc/min 
V 2 rr 2.005 x 10 4 cc/min 


(T = 2,073°K, p = 0.5 atm) 


( 49 ) 

(50) 


and the corresponding Reynolds numbers were calculated to be 

% ei = 56.2; N Be2 = 37 (51) 

Under these new flow conditions Vp for the central jet was calculated 
as v p = 179.5 cm/sec. Since p = v s /v p and downstream interaction 
between central and peripheral jets occurs, it is obvious that p 
depends on x. An average value of p was estimated on the basis of two 
limiting values of v s = 17.6 cm/sec and v s = 134 cm/sec; for pAvg =0.223 
the total velocity at the centerline, (v x ) 0 = v"i + v s , was calculated as 


( v x)o 


= 1,180 + 40 

x 


(52) 


and the jet width as 


W = 0.24 (x)°- 777 


(53) 


where (v x )o [=] cm/sec; x, w, [=] cm. The results are shown in Fig. 19 . 
Similar calculations for the peripheral jets and available data, 34 gave 
x c = 8.5 cm, p = 0.6l and 

(-x)o = aa. + 82 (54) 
W = 0.294 (x)°- 39 . (55) 
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X, Cm 

Fig. 19 - Velocity at the Center Line and Width Variation for Downstream 
of Central Jet 


2W Cm 


The results are shown in Fig. 20. The rms value u' was calculated using 
available data. 35 For the central jet and at a radial distance of 
r = 1 cm, 



u , = 1.8 x 10 2 
(x - 1.02) 



[=] cm/sec 


(56) 


where x[=]cm, and the function g was obtained from reference [ 36 ]. 

Since those data were for much higher jet Reynolds number, the predicted 
values of u' from equation (52) are rather high. Due to lack of other 
data this calculation was retained for comparison. Hinze 37 gives the 
following correlation for the axial (x-direetion) eddy viscosity, (e^xx. 

( e m)xx = 0.013 v p d; [=] cm 2 /see (57) 

where Vp [=] cm/sec and d [=] cm. 

The corresponding Reynolds (eddy) stresses are given by the 
well known expression 

^x^Sc = " ^ G m^xx ^ (58) 


Using Figs. 19 and 20 or equations (48) and (50), the velocity gradients 
at the centerline may be calculated.. The following expressions for u* 
for both the central and peripheral jets can be obtained. 


Central jet: u ! = 52 

(59) 

X 

Peripheral jets: u' - 20. 3 

( 60 ) 


x 


where x[=]cm, u'[=]cm/sec. These results are summarized in Fig. 21. 

It is worth noticing that u' for the peripheral jets is lower than for 
the central jet. Half-way in the x-direction of a ~35 cm long reaction 
zone, u' - 1.3 cm/sec; the velocity (v x ) 0 at the same x is 100 cm/sec 
[Fig. 20] or u' - 0.013 (v x ) Q . 
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Fig. 20 - Peripheral Jets 




From Fig. 21, one may graphically calculate a reactor length 
for a chosen intensity of segregation. Additionally, an analytic solu- 
tion was attempted. For gases, with the Schmidt number very close to 
unity, equations (39) and (40) give for the time constant, x, 



for r c = 1 inch =2.54 cm, and for the central jet (equation 24) 
x = 0.022 x, x [=] cm and x [=] sec. The intensity of segregation 
becomes 


I s (x,t) = exp 


t 

0.022x 


( 62 ) 


The total volumetric flow rate (Streams 1 and 2) was 
V = 1.373 x 10 4 + 2.005 X 10 4 = 3-378 X 10 4 cc/min = 5.62 x 10 2 cc/sec 

(63) 


From the geometry shown in Fig. 22, it is easily seen that R = 2.54 
+ (l.27/L)x [=] cm. If v = v(x) is the average velocity based on the 
total cross sectional area, 

_ dx V [=] cc/sec _ 5§2 

dt rtR2 ^2.54 + h21 xj 2 


If L = 35 cm, equation (64) represents the intensity of segregation as 
I s (x) = exp[-0.279 (4.4 • 10 -3 x 2 + 9.2 • 10" 2 X + 6.45)] = e -4 

(65) 

where x [=] cm and 8.9 < x < 35. Table XXI lists the values of I s (x) 
from equation ( 65 ). 
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Fig. 21 - Variation of Axial 
With X 






Fig. 22 - Geometry for Calculating Residence Time 
in Reactor 


Table XXI - Intensity of Segrega- 
tion from Equation ( 65 ) 


x[=]cm 

t 

Is = e-Y 

9 

2.13 

0.119 

10 

2.16 

0.115 

15 

2.47 

0.094 

20 

2.78 

0.072 

25 

3.21 

o.o4o 

30 

3.77 

0.023 

35 

4.20 

0.016 
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It is obvious from Table XXI, that the calculation of time based on the 
total cross-sectional area gives high values for t, thus decreasing the 
intensity of segregation very rapidly. 

Recent research 38 in gas jet mixing suggested that the length 
of the potential core x c - L s - r 0 , at least for the initial period of 
mixing i.e., for 0 < x < x c . A conservative calculation performed on 
this basis resulted in the following. For the central jet and 
0 < x < x c = 8.9 cm, r c « x c , u' - 7 cm/sec, the u' close to the end 
of the potential core [Fig. 21], the time constant, t, was calculated. 

t = (0.341) ( 8.9 cm ) = o. 434 sec (66) 

7 cm/sec 

For the region: 8.9 < x < 35 cm, r 0 =* 2.5 inches = average of the two 

diameters of 2" and 3". u' =* 2 cm/sec (i.e., u' at the end of the 
reaction zone) 

t = (0.341) (6.35 on) _ 1.07 sec ( 67 ) 

2 cm/sec 

The results for the intensity of segregation and time are shown in 
Fig. 23. A time for the central jet to travel the total length of the 
reactor was calculated. 

t = 3 * * * * 8 -2. c . m . + fe j eg = 0.33 sec (68) 

179 cm/sec 93 cm/sec 
(see Fig. 19) 


The corresponding intensity of segregation from Fig. 23 was I s s 0.66, 
or 34% mixing. 


3. Reactor Mixing for Zn 2 Ti0 4 

The proposed CaW0 4 reactor mixing performance is estimated for 
the production on Zn 2 Ti0 4 according to the reaction: 

2ZnCl 2 + TiCl 4 + 4H 2 + 4C0 2 = Zn 2 Ti0 4 + 8HC1 + 4C0 ( 69 ) 

The reactor configuration, gas flow rates, operating temperature and 
pressure are given in Fig. 24. The amounts of 0 2 and H 2 in stream (l) 

are the ones for the oxy-hydrogen flame and for the purpose of this 

estimation they were based on heat supply of 40 kcal/hr by the oxy- 
hydrogen flame. 
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Fig. 23 - Mixing of Central Jet 


6l 
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Gas flow rates 


Stream (l): 


ZnCl 2 

4oo 

cc/min 

(STP) 

h 2 

248 

t! 

It 

C0 2 

4oo 

ft 

It 

0 2 

124 

ft 

It 

ream (2): 

fiCl 4 

200 

cc/min 

(STP) 

h 2 

800 

ft 

tt 


Temperature: 1,300°C 

Pressure: 0.5 atm. 

The volumetric flow rates of streams (l) and (2) at 1,300°C 
and 0.5 cc/min were calculated as: 


V L = (1,172) x (^ff ) - 129480 cc / mln (70 ) 

V 2 = (1,000) x ( -i- ) ( 1^598 \ s 10,700 cc/min (7l) 

'0.5 ' ' 298 > 

The corresponding velocities of the central jet, Vi, and each of the 
eight peripheral jets were calculated (at the operating temperature and 
pressure) 


( 12,480 cc/min) 
(2.53 cm 2 ) (60 sec/min) 


82.3 cm/sec 


(72) 


v 2 (8,460) _ cm/sec 

(0.315) (60) (8) 


(73) 


The kinematic viscosities of streams (l) and (2) were estimated as 
01 = 1.77 cm 2 /sec and o 2 = 1-38 cm 2 /sec (see Table XX) respectively. 
The corresponding Reynolds number of the central jet was calculated as 

N Re = 59- 
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Jet mixing estimation 


Equation (45) with a ratio of velocities q = v s /vp 
= 56/82 = 0.68, gives for the length x c of the central jet core region. 

x c = (4 + 12 x 0.68) x (1.27) = 15.4 cm (74) 

The velocity at the centerline of the jet (v x ) 0 was calculated from 
equation (46) as: 


( v x)o “ Vs x 


v p - v s 


_c 

X 


or 


( v x)o = ~ + 56; x[=]cm, (v x ) 0 [=]cm/sec 


(75) 


Equation (57) for the eddy viscosity (ej^)^ gave: 

(em)xx = (0.013) (82) (1.27) s 1.36 cm 2 /sec 
From equation (75), 

d ^ v x^o _ _ 400 
dx x 

which combined with equation (76) reduced equation .(52) to: 


(76) 


(77) 


u ' = yj + (1-36) [ = ]cm/sec; x[=] 

X^ X 


cm 


(78) 


The time constant of mixing, t, was calculated from equations (39) and 
(40) at a radial distance r Q = 1 in.. 


T 


H 


3 x lOjUlj (2.54)T 0>93 sec 


1.4 


(79) 
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The intensity of segregation 

I s = e ”V°- 93 , t [=] sec (80) 

is plotted vs mixing time in Fig. 25- At a downstream distance x = 20 cm, 
the velocity at the centerline of the central jet is: (v x ) D •= 76 cm/sec 

(from equation 75)- This velocity was used to calculate a residence 
time in the reactor. This is a conservative estimate as far as mixing 
because the jet velocity drops continuously downstream, and the central 
jet is the faster one; in addition the velocity profile has its maximum 
at the center line of the jet reaching the value v s = 56 cm/sec at the 
most, towards the boundaries of the expanding nozzle of the reactor. 

A mixing time based on the Lm 35 cm reactor length and the volume 
velocity was calculated as t = (35 cm)/(76 cm/sec) - 0.46 sec, giving 
I s - 0.60 for the intensity of segregation, i.e., 40% mixing. Hence, 
even this conservation estimate gives an adequate degree of reactant 
mixing in the case of Zn 2 Ti04 production. 
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III. EXPERIMENTAL PROCEDURE 


Serious experimental limitations of previous investigations were 
reviewed and analyzed. Methods for alleviating or circumventing these 
problems were based on improved component and process designs. The 
results of theoretical gas mixing studies were incorporated, as were 
new devices such as induction heated reactors and sonic agglomerators . 


1. Reactor 

A silicon carbide tube furnace 21 inches long and 2.25 inches 
ID was used to heat a high-purity, recrystallized A1 2 0 3 tube 36 inches 
long with 38 mm ID and 46 mm OD. Reaction of gases and subsequent 
powder production occurred in a 6 inch hot zone. Use of a silicon 
carbide heating element provided a maximum temperature capability of 
1500°C . Temperature was measured with a Pt-Pt 10$ Rh thermocouple 
positioned between the A1 2 0 3 tube and the SiC heating element in the 
center of the hot zone. The furnace was positioned vertically with the 
A1 2 0 3 reaction tube passing through the center and extending four inches 
above the heating element . The tube furnace used to heat the zinc and 
zinc oxide extended below the SiC heating element. 

Gases were introduced into the reaction zone at the bottom and 
passed up through the furnace. When zinc was used as the zinc vapor 
species, all other gases were introduced through a A1 2 0 3 injector nozzle 
(6 mm xlO mm) central to the zinc crucible and extending within one inch 
of the hot zone. The A1 2 0 3 reaction tube was sealed with a stainless 
steel plate and a silicone rubber gasket. The top of the A1 2 0 3 reaction 
tube (outlet) was connected to a stainless steel tee joint and sealed 
with a silicone rubber gasket. The branch of the tee was connected to 
a rubber vacuum hose and the run of the tee was used as a sight window. 

Three reaction furnaces were designed and constructed to extend 
the range of system conditions for powder production. Intermediate 
temperatures to 1500°C were provided by a resistively-heated silicon 
carbide furnace. The silicon carbide element was thermally insulated 
with K-30 insulating fire brick. One-half of the symmetrical temperature 
profile, with maximums of 1050°C and 1200°C is shown in Fig. 26. Taking 
the effective isothermal hot zone as that length within 50°C of the 
maximum, a hot zone of approximately ten inches (25 cm) resulted. 

Temperature capability was extended beyond l800°C with a 
molybdenum-wound, hydrogen-protected furnace, Fig. 27. The 60-mil 
molybdenum wire was coiled about a wax core. Successive layers of 
Alundum cement containing 20 volume per cent of polyvinyl alcohol were 
built up on the molybdenum element, to a total thickness of 1-1/2 inches. 
The wax was melted, leaving a hollow tube of Alundum cement with the 
molybdenum wire firmly embedded in it. This element was centered in an 
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Temperature (°C) 



Distance from Furnace Center (inch) 

Fig. 26 - Temperature Profile of Silicon Carbide Resistance- 
heated Furnace 
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Fig. 27 - 22-inch Molybdenum-Wound Furnace 



18-inch diameter drum, and the remaining volume filled with bubble 
al umi na insulation. Electrical leadthroughs in the outer shell served 
as connections to the power supply. The drum was sealed and checked 
for gas -tightness. 

The molybdenum element was protected from oxidation during 
firing by hydrogen. The polyvinyl alcohol burned out during the first 
firing, providing porosity in the core for flow of hydrogen over the 
element. Power was supplied by a variable transformer, and temperature 
was monitored with an optical pyrometer . Furnaces with four and seven 
inch isothermal zones were constructed. 

Existing reactors heated by silicon carbide and molybdenum 
elements were modified to the configurations shown in Figs. 28 and 29. 
Respective temperature profiles of these reactors are shown in Figs. 26 
and 30. In both instances, the halide is produced by passing chlorine 
through a bed of appropriate oxide. Figure 31 shows the arrangement of 
molybdenum-heated reactors with pumps and flow controls . 

Operating temperature was extended beyond l800°C by the con- 
struction of a high frequency, induction -heated reactor. Fig. 32. The 
temperature profile of the reactor is shown in Fig. 33- The arrangement 
of flow controls, pumps, furnace and sonic agglomerator are shown in 
Fig . 34 . 


The reactor contained an internal chlorinator for zinc oxide 
or zinc orthotitanate opening into a 13 -inch long diverging nozzle which 
was the reaction zone. The chlorinator and nozzle were separated by an 
alumina disc. The reactor chamber was a three -inch diameter alumi n a 
tube. Radiation losses were minimized by two concentric molybdenum 
heat shields, .010" each. The reactor and heat shields were positioned 
within a 6. 5- inch diameter fused silica tube. Argon flows prevented 
heat shield oxidation. 

The induction power supply was a 15 kW, 10,000 Hz, motor 
generator feeding an eight tap autotransformer in a portable work 
station. The induction coil was 7 inches in diameter, 20 inches long, 
and contained 18 evenly spaced turns of 5/8" diameter copper tubing. 


2. Powder -Gas Separation 

A survey of the literature identified several methods of solid- 
gas separation. 39 " 50 Four separation processes were considered; namely, 

(1) centrifugation, 

(2) electrostatic precipitation, 

(3) thermal precipitator, and 

(4) cyclone. 
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Fig. 29 - Modified Molybdenum-Wound Furnace 





Fig. 31 - Molybdenum-Heated Reactors with Pump and Flow Controls 
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Fig- 34 - Induction System 




The performance of a particle collector is described by col- 
lection efficiency. Regardless of the figure of merit, performance is 
not a specific characteristic of a given collector but depends on the 
physical properties of the dispersoid. Figure 35 shows characteristics 
of particles and dispersoids together with the applicable types of gas 
cleaning. The primary distinguishing feature of gas dispersoids is 
particle size. The most widely used unit of particle size is the 
micron (p). The particle size of a gas dispersoid is usu ally taken as 
the average or equivalent diameter of the particle. 

The most widely used type of dust -collection equipment is the 
cyclone, in -which dust-laden gas enters a cylindrical or conical chamber 
tangentially at one or more points and leaves through a central opening. 
The dust particles, by virtue of their inertia, w ill move toward the 
outside separator wall from which they are led into a receiver. At 
operating conditions commonly employed, the centrifugal separating force 
or acceleration may range from 5G, five times gravity, in very large 
diameter, low-resistance cyclones, to 2500G in very small, high-resistance 
units . In a cyclone, the gas path involves a double vortex with the gas 
spiralling downward at the outside and upward at the center. When the 
gas enters the cyclone, its velocity undergoes a redistribution so that 
the tangential component of velocity increases with decreasing radius 
as expressed by V c t ~ r“ n . A cyclone is essentially a settling chamber 
in which gravitational acceleration is replaced by centrifugal acceler- 
ation. 5 1 


Various centrifugal collectors or cyclones have efficiencies of 
only 8-40$ at one micron and fall rapidly with decreasing particle size. 
Electrostatic precipitation is limited by the high resistivity of the 
powder and the increased arc-over at reduced operating pressures. Of 
the remaining techniques, sonic agglomeration offers the greatest promise 
for laboratory collection. Construction is simple, even for corrosive 
service, and collection efficiency is independent of temperature. 

A sonic agglomerator subjects the particle -laden gas stream to 
a standing accoustic wave. Vibration causes particle migration to the 
antinodes where flocculation occurs . Sonic flocculation is thought to 
result by a combination of: (a) co-vibration of the particles and gas, 

(b) sonic radiation pressure and (c) hydrodynamic forces between 
particles . The degree of flocculation is predictable from the geometry 
of the tube, the power input and the residence time of the gas stream. 

The floes should be of sufficient size to permit collection by conven- 
tional means, as in a cyclone. 

To determine the geometry of an agglomeration tower, an index 
of agglomeration, I a , was defined such that 

la = f ( 81 ) 

o 

where 
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Fig- 35 - Characteristics of Particles and Particle Dispersoids. 

(Courtesy of the Stanford Research Institute, prepared 
by C. E. Lapple ) 51 
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d = final agglomerate diameter 
d D - initial particle diameter 

An index of 100 was necessary to change 0.1 micron powder into 
10 micron agglomerates which can be collected by a cyclone. From Boucher 52 
an index of agglomeration of 100 required a residence time in the tower 
of about ten seconds. For typical flow rates of 100 cc per second, a 
tower volume of 1000 cc was needed to satisfy this requirement. Thus, 

1000 = 5f2i ( 82 ) 


where 


D = tower diameter 
L = tower length 

A diameter of 5-1 cm (2 inches) and length of 42 cm (l6.5 inches) were 
chosen as the tower dimensions. 

To determine the sonic power requirements, the coalescence 
constant, C, was determined. It is related to the index of agglom- 
eration by 


I a = e<«/3> (83) 

where t is the residence time. For I a = 100 and t = 10 seconds, the 
coalescence constant, C, was 1.4. Boucher has shown that the average 
sound intensity, I, within the tower is: 

I = C 2 watts/cm 2 (84) 

The cross-sectional area of the tower, A, is 20 cm 2 . Thus, 


Power (P) = IA = 40 watts 


(85) 


Based on the design considerations a sonic agglomerator , Fig. 
36, was constructed. Sonic input was supplied to 11,000 Hz and power 
levels to 60 watts . A flexible ^lar film protected the generator from 
the gas stream and kept the system closed. A glass cyclone, Fig. 37 was 
placed downstream to collect the agglomerates . 

Gases exiting the cyclone were passed through a liquid nitrogen 
cold trap to remove condensable vapors and prevent excessive p um p oil 


80 







19mm 



10 cm 




contamination. Non-condensable vapors were vented to a hood. For con- 
tinuous operation product gases would be separated, probably with a wet 
scrubber, and recirculated to complete the loop. 


3* Gas Scrubbers 

An on-line vapor phase reaction system should be equiped to 
separate, recover, or convert into useful products the gaseous HC1, Cl 2 , 
and CO. A review of the literature to identify efficient and practical 
means of recovering product gases indicated some of the wet-scrubbing 
techniques to be 33_6? 

spray towers, 
deflection washers, 

mechanical (combination of dry and wet), 
atomizing (venturi throat-type), and 
wetted fibers and packed towers (only method 
described for acid mist removal). 

in which a liquid is employed to achieve or assist in the removal of 
powders from gases. Figure 38 shows a block-type flow plan and a 
schematic of a typical laboratory-type scrubber. The basic and most 
essential part of this scrubber is the absorption tower, consisting of 
water flowing over activated charcoal. 

The flows of the two fluids (water down and gas up) in the 
absorption tower are important. The water must flow at a minimum wetting 
rate (MWR) so that it wets all the charcoal. The relative flow bf the 
two fluids (i.e., water and gas) must not be so large, however, that 
the upward flow of gas impedes the downward flow of water by bubbling 
through it, a point known as the flooding rate is reached, at which the 
flow of water is impeded. The tower is designed to operate at 0.2°C 
and remove HC1 and Cl 2 from the gas stream. The HC1 is removed by 
direct reaction with water. The absorption coefficient (that is, the 
solubility of HC1 in water at 0°C and at l/30 atm.) is 2.74 g/100 ml. 

The Cl 2 is removed by the reaction 63 

2 Cl 2 + 2 H 2 0 + C 4 HC1 + C0 2 . (86) 

Because the solubility of Cl 2 in H 2 0 is small, (i.e., ~ 0.1 
g/100 ml at 0°C and l/30 atm.) the above reaction will be the dominant 
mechanism for the removal of Cl 2 from the gas stream. With the HC1 and 
Cl 2 removed for the gas stream, the gas out of the absorption tower 
consists of water vapor, CO, (from the original product gas) and C0 2 
(from the Cl 2 reactions). At this point the gas stream is dried (pos- 
sibly by 2 CaS0 4 drying column), and the C0-C0 2 mixture is passed through 
a separator. 
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Using available data for gas imputs in the vapor phase-growth 
apparatus, approximately 4.3 x lO 5 * 1 molecules/minute of solid species 
are produced. For example, for Zn 2 Ti 04 this corresponds to 4 x 4.30 x 
lO 21 molecules of HC1. By converting the number of HC1 molecules to 
moles and then to grams, we get approximately 1.0 grams. If the solu- 
bility of HC1 in water at l/30 atm and 0°C is 2.74 gm HCl/100 ml water, 
then 1.04 gm dissolves in 38 ml of water. If an equal amount of HC1 
is produced by the Cl 2 reaction in the absorption tower, approximately 
80 ml/min of water is necessary to remove the HC1 and Cl 2 from the gas 
flow. Using a conservative safety factor (and realizing that the flow 
rates assumed in this treatment are high) 200 ml/min of water should 
be sufficient. The cold trap and back-flow trap upstream from the 
absorption tower are added variations that prevent solid particles from 
entering the absorption tower. 


4 . Chlorinator 

A chlorinator for titanium was designed, fabricated, and 
installed on the apparatus . The chlorinator was electrically insulated 
by wrapping pressed Fiberfrax around the surface; Kanthal wire was 
wrapped over the Fiberfrax. The chlorinator was thermally insulated 
by wrapping the Kanthal wire-wrapped chlorinator with Kaowool blanket- 
type insulation. The chlorinator temperature profile was stabilized 
and calibrated. Because of a temperature drop at the outlet end of the 
chlorinator, the Kanthal wrapping was modified. The rewrapped Kanthal 
was stabilized and calibrated. 

A double-stage chlorinator. Figs. 39 and 40, was constructed 
for the production of Zn 2 TiC> 4 . The purpose of the double-stage chlorin- 
ator was to chlorinate zinc oxide and titanium dioxide, and produce a 
homogeneous mixture of zinc chloride and titanium tetrachloride before 
the chlorides entered the reaction zone . 

To establish a controlled halide flow to the reaction chamber, 
the conversion of chlorine to metal chloride must be determined. To 
measure the conversion efficiency, a chlorinator was constructed as 
shown in Figure 4l. The Inconel chamber containing the metal was fitted 
with thermowells for accurate mapping of radial and longitudinal temper- 
atures. External heating with asbestos -insulated Nichrome wire main- 
tained the chlorinator and transfer lines at the desired temperature. 

A Pyrex condenser immersed in liquid nitrogen was used to collect the 
chloride. The conversion efficiency was calculated as the ratio of the 
actual weight of chloride in the condenser to the theoretical weight 
for 100fo conversion. 
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Fig. 39 - Double-Stage Chlorinator 
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Fig. 40 - Double -Stage Chlorinator 
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5. Gas Temperature Measurement 

Because prediction of reaction products using thermodynamic 
data is based on constant temperature, it was necessary to know the 
actual temperature of the gases to verify reaction temperature. An 
alumina tube (10 mm in diameter) was positioned in the center of the hot 
zone with a thermocouple placed next to it. The temperature of the alum- 
ina tube was obtained by sighting an optical pyrometer on a mirror reflec 
tion through the viewing port while at the same time recording the thermo 
couple temperature. This provided a calibration of the radiation loss 
due to the mirror and glass viewing port. The thermocouple was next 
positioned between the alumina reaction tube and the silicon carbide 
heating element at the center of the hot zone at the same height as the 
10 mm alumin a sight tube. The system was evacuated and temperature 
measurement of the interior alumina sight tube and the thermocouple out- 
side the reaction chamber were recorded for a range of temperatures, 
pressures, and gas flow rates. Hydrogen, carbon monoxide, and carbon 
dioxide were used as flow gases in the temperature measurements because 
they would not form solid reaction products which would interfere with 
the optical pyrometer readings. The effects of gas flows and pressure 
on the reaction zone temperatures are shown in Fig. 42. The interior 
of the reaction zone was 10°C lower than the exterior when pressure was 
at a mini mum and no gases were flowing. The interior temperature 
decreased as gas flow rates were increased. At the highest flow rate 
the interior temperature was approximately 40°C lower than the exterior 
temperature range of 1000°C to l400°C . A temperature profile of the 
vapor reaction furnace is shown in Fig. 43. 


6. Gas Flow Control 

Reactant gas flows were monitored using low- flow rotameters . 
Chemically pure CO, welding grade C0 2 , and 99-97 per cent H 2 were passed 
through calcium sulfate drying columns to remove moisture before moni- 
toring. Chlorine (99-965 pure) was introduced into the system at three 
positions through three separate flowmeters and transfer lines. Refer 
to Figs. 44 through 46 for vapor phase apparatus, gas control panel, and 
a schematic of the system gas flow control, respectively. 

To provide greater flexibility in gas flow rate ranges, two 
sets of flow meters were operated at different pressures. The low- flow 
rotameters were calibrated in cc/min at 0 psig and 5 psig for H 2 , CO, 
and C0 2 , and in cc/min at 0 psig and 25., psig for Cl 2 . For calculation 
of desired flow rates, gas pressure through the rotameters had to be 
considered. Gas pressures through the rotameters and flow rates were 
adjusted by micrometer needle valves below the rotameters. 

System pressure was controllable from atmospheric pressure to 
5u using two Duoseal vacuum pumps in parallel and a bleed valve which 
opened to the atmosphere. Pressure was monitored with two absolute 
diaphragm gages which had ranges of 0-50 mm Hg and 0-760 mm Hg. 
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Fig. 42 - Actual Temperature Inside Reaction Tube vs Temperature 
Outside Reaction Tube for Various Gas Flow Rates 
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Fig. 43 - Temperature Profile of Vapor Phase Reaction Furnace 
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Fig. 44 - Vapor Phase Apparatus 
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Fig. 45 - Gas Flow Control Panel 
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Fig. 46 - Vapor Phase System 














7. Operation Procedure for Vapor Phase Runs 

The source metals or metal oxides were weighed and placed in 
the reaction zone and chlorinator. The system was evacuated and all 
portions of the system were brought up to temperature and stabilized 
before each run. Liquid nitrogen was poured into the cold trap. Gas 
flows were initiated in the following order: Cl 2 to metal or oxides, 
C0 2 , and H 2 . At the completion of a run gases were stopped in the 
reverse order from initiation, the system was closed above the cold 
trap, and the cold trap was removed and placed in a hood while the 
trapped gases evaporated. Once the chlorinator cooled to room tempera- 
ture, the reaction tube was removed and the reaction products analyzed. 
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IV. EXPERIMENTAL RESULTS 


A. ALUMINUM OXIDE (AI2O3) PRODUCTION AND CHARACTERIZATION 
1. Aluminum Oxide Production 

Aluminum oxide powders were prepared by the following reaction: 

2 AICI3 + 3 H a + 3 C0 2 ^ A1 2 0 3 + 6 HC1 + 3 CO (87) 

The reaction free energy and equilibrium constant were cal- 
culated over the temperature range 500-2000°K. The thermodynamic data 
are summarized in Table XXII . 

The growth conditions for the formation of aluminum oxide are 
given in Table XXIII. Aluminum chloride was generated as discussed 
previously. Weight-loss measurements on the source aluminum indicated 
a minim um chlorination efficiency of 80-95 per cent for the various runs. 

Reactant gases were admitted to the reaction chamber through 
a central injector nozzle. The reaction chamber was a horizontal 
alumina tube, 38 millimeters in diameter. The tube was heated in a 
molybdenum-wound, hydrogen-protected furnace . Gases were reacted in a 
10-cm long isothermal zone. Aluminum oxide powders nucleated and grew 
during residence in the isothermal zone, beyond which they deposited 
on the tube walls. Samples were carefully selected from various regions 
within the tube for characterization. 


2. Aluminum Oxide Characterization 

The crystalline phases in the powder samples were identified 
with x-ray diffraction (see Appendix D). In all samples, the principal 
phases were transition aluminas, predominantly gamma and delta. Alpha 
alumina (corundum) was present as a secondary phase. Traces of 
unreacted aluminum chloride were observed in a few samples collected 
near the cold end of the tube where the temperature would permit 
deposition. 

The particle size of the powder samples was determined by 
lineal analysis of electron micrographs. Several typical micrographs 
are shown in Figs. 47-50. 

Comparing the particle size with the growth conditions for each 
run, several trends appear. 

In general, increasing system pressure decreases the gas velocity 
in the reaction chamber, allowing first a longer time for gas mixing and 
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Table XXII - Free -Energy Changes and Equilibrium Constants as a Function 
of Temperature for Aluminum Oxide Production 


T 

( # c) 

T 

(°K) 

AG° 

Reaction* 

(kcal/mole) 

K p 

Reaction* 

227 

500 

-63.3 

3.3 x 10^ 

727 

1000 

-75*8 

3.0 x 10 16 

927 

1200 

-80.2 

3.4 x 10 14 

1127 

1400 

-85.4 

1.8 x 10 13 

1227 

1500 

-87.6 

5.0 x 10 12 

1327 

1600 

-90.5 

2.0 x 10 12 

1427 

1700 

-92.6 

6.9 x 10 11 

1527 

1800 

-94.6 

2.7 x 10 11 

1627 

1900 

-97.4 

1.4 x 10 11 

1727 

2000 

-100.0 

7.4 x IO 10 

*2 AICI3 + 3H2 

+ 3 CO ^ 

A1 2 0 3 + 6 HC1 + 3 CO 
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Table XXIII - Conditions for Aluminum Oxide Powder Production 


Flow Conditions (cc/min) p Total Temp. Particle 

1 *— — ■ — ■ 1 ■ i— ..i . 4m ■" ■ i . » . » / \ 


Run No. 

h 2 

co 2 

CO 

C12XS 

Cl 2 A 1 

(Torr) 

(°c) 

Size (p) 

VA -2 

1000 

200 

- 

- 

30 

15 

1500 

~ 0.01 

VA -3 

1000 

200 

- 

- 

30 

200 

1500 

~ 0.01 

VA -4 

1000 

200 

- 

- 

4o 

350 

1500 . 

~ 0.01 

VA-6 

950 

800 

100 

20 

40 

100 

1750 

0.05 

VA-8 

950 

8oo 

100 

20 

70 

340 

1700 ) 









1750 ( 

1800 ( 

0.07 








1850 ; 


VA -9 

950 

800 

100 

20 

70 

260 

1850 

0.18 

VA -10 

6oo 

500 

210 

45 

5 ° 

660 

1800 

0.13 

VA -13 

100 

100 

- 

- 

80 

50 

1700 

0.02 

VA -15 

100 

100 

- 

- 

80 

100 

1750 

0.09 

VA -16 

50 

50 

- 

- 

40 

50 

1750 

0.05 

VA -17 

100 

100 

- 

- 

80 

200 

1750 

0.12 

VA- 1 B 

50 

50 

- 

- 

40 

50 

1600 

0.14 



Fig. 47 - Electron Micrograph of Alumina Powder, 67,200X 
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Fig. 50 - Electron Photomicrograph of Alumina Powder, 43,4-OOX 
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heat transfer in the reactor entrance, and second a longer residence 
time for powder growth in the hot zone. Increased heat transfer to the 
gas stream and increased mixing should result in increased efficiency 
of conversion. This was substantiated by the increased amount of powder 
produced at higher pressures. In addition, the increase in residence 
time in the hot zone should result in an increase in average particle 
size. The increase in size from 0.09p for Run No. VA-15 to 0.12p for 
Run No. VA-17 verified this assumption. 

To determine the effect of temperature upon particle size, 
temperature was varied in VA-8 from 1700°C to l850°C, in 50°C incre- 
ments . The particle size distribution of the resulting powders was 
much wider than that obtained at fixed conditions. A range from 0.01 
to 0.20 was observed with an average particle size of 0.07q. From a 
probability plot of particle size versus per cent finer (Fig. 51 ), the 
particle size distribution appeared to be a sum of several individual 
normal distributions. The results of VA-9 suggest that the large size 
of VA-8 can be attributed to the higher temperatures. However, com- 
paring VA-16 with VA- 18 , it is apparent that a decrease in temperature 
resulted in considerable increase in particle size for these particular 
input gas conditions . The effect of temperature on particle size is 
not precisely defined from the results of these experiments. 

Comparing the results of VA-15 with VA-16, the effect on 
particle size of the number of moles input is apparent. Doubling the 
inputs and simultaneously doubling system pressure to maintain the same 
gas velocity results in approximately a two- fold increase in particle 
size from 0.05 to 0.09|i. Thus growth appears to be directly related to 
the availability or concentration of the growth species . 


B. TITANIUM DIOXIDE (Ti0 2 ) POWDER PRODUCTION AND 
CHARACTERIZATION 

1. Titanium Dioxide Powder Production 

The powder was formed by the following chemical reaction and 
the subsequent homogeneous nucleation of rutile: 

TiCl 4g + 2 H 2g + 2 C0 2g * Ti0 2g + 4 HCl g + 2 C0 g . (88) 

The reaction free energy and equilibrium constant for the 
reaction (Table XXIV) show thermodynamic feasibility over the entire 
temperature range. Titanium chloride was generated by controlled 
chlorination of titanium sponge at 400°C, and transferred to the reaction 
chamber at 300°C to prevent condensation in the lines. Runs were made 
over the temperature range 800-l400°C. The silicon carbide furnace 
heated the reaction chamber. Gases were injected into the chamber 
through a five millimeter (inside diameter) alumina tube which extended 
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Fig. 51 - Particle Size Probability Plot of Alumina Powders 
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Table XXIV - Free-Energy Changes and Equilibrium Constants 
as a Function of Temperature for Titanium 
Dioxide Production 


T 

(°c) 

T 

(°K) 

AG® 

Reaction * 
(kcal/mole)' 

Kp 

Reaction* 

27 

300 

- 4.60 

2.25 x 10 3 

127 

4oo 

-10.55 

5.82 x 1CP 

227 

500 

-17.15 

3.14 x 1(7 

327 

600 

-18.20 

4.27 x 1CP 

427 

700 

-21.85 

6.64 x 1C P 

527 

800 

-25.45 

8.98 x 1CP 

627 

900 

-28.90 

1.04 x 10 7 

727 

1000 

-32.10 

1.04 x 10 7 

827 

1100 

-35-55 

1.16 x 10 7 

927 

1200 

-39.39 

1.49 x 10 7 

1027 

1300 

-42.30 

1.29 x 1CF 

1127 

1400 

-45.90 

1.46 x 10 7 

1227 

1500 

-49.15 

1.45 x 10 7 

1327 

1600 

-53.25 

1.88 x 10 7 

1427 

1700 

-56.45 

1.81 x 10 7 

1527 

1800 

-60.60 

2.28 x 10 7 

1627 

1900 

-64.50 

2.63 x 10 7 

1727 

2000 

-67.40 

2.32 x 10 7 

*TiCl 4 + 2 H2 

+ 2 C0 2 Ti0 2 

+ 4 HCl + 2 CO 
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within two inches of the hot zone.. Residence time was controlled with 
system pressure and was calculated for a 25 -centimeter isothermal zone. 
System conditions and particle sizes for specific runs are listed in 
Table XXV. Analysis by x-ray diffraction of titania powders showed 
only rutile present at 1200°C and l400°C . A minor amount of anatase 
was discovered in the 1000°C powders, and at 800°C specimens were about 
equal portions of rutile and anatase. The presence of anatase was con- 
firmed by the appearance of the typical anatase morphology in powder 
micrographs, Figure 52. The dominance of pyramidal faces over prisms 
in anatase makes its octahedral habit readily distinguishable from the 
prismatic rutile modification. 


2. Titanium Dioxide Powder Characterization 

X-ray diffraction analyses of the powder from both runs, using 
both diffractometer and powder camera methods, showed the material to 
be rutile. No other phases were detectable. Typical micrographs of 
rutile powders. Figure 53-56, describe the range and uniformity of size 
obtained. Average particle size ranged from 0.06 to O .38 micron. 
Distributions were narrow, Figure 57, comparing favorably with or 
superior to materials from other processes. Considerable improvement 
in surface roughness over commercial oxidation process rutile was 
observed. Surface areas from 2.3 to 6.6 square meters per gram resulted 
in surface roughness values of 1.3 to 1.7* This improved surface state 
should enhance the stability of the pigment in extraterrestrial thermal 
control coatings . 

The effect of residence time and hot zone temperature on the 
particle size of vapor-grown rutile is presented in Figure 58 . For 
three-dimensional growth, particle volume increased with residence time 
in the isothermal zone. Since the cube of the diameter provides some 
measure of particle volume, this parameter was chosen for the abscissa 
in Figure 58 . It was apparent that particle volume increased linearly 
with residence time, indicating a constant volume growth rate for a 
constant temperature and reactant gas input. The variation in particle 
size for very short residence times was indeterminate, but a finite 
ordinate intercept would be anticipated, representing an induction time 
required for reaction and nucleation of the growth phase. The volume 
growth rates, or slopes, from Figure 57 were strongly influenced by 
temperature, Figure 59. A maximum in growth rate at 1000°C may be 
explained by considering separately the effects of nucleation and 
growth. ' Nucleation rate was qualitatively observed to decrease with 
temperature. Profuse powder production at l400°C decreased slighly at 
1200°C and was dim inished still further at 1000°C . Very little powder 
was produced at 800°C. High nucleation rate may have depleted the gas 
stream of the growth species, and little material remained for growth 
at the higher temperatures. As reaction and nucleation rates decreased 
with temperature, fewer nuclei were produced and their resultant growth 
was less inhibited by concentration depletion. The net result was a 
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Table XXV - System Conditions for Ti0 2 Pigment Production 


Run No. 

Temp . 

C°c) 

Pressure 
(mm. Hg) 

Gas Velocity 
(cm/sec) 

Residence Time 
(sec) 

Average Particle 
Size (microns) 

T- 6 

1400 

50 

62.2 

0.4 

0.09 

T- 4 

1400 

200 

15.5 

1.6 

0.18 

T- 5 

1400 

400 

7.8 

3.2 

0.18 

T- 7 

1400 

660 

4.7 

5.3 

0.16 

T-18 

1400 

760 

4.1 

6.1 

0.19 

T-12 

1200 

50 

54.8 

0.45 

0.06 

T- 8 

1200 

100 

27.4 

0.91 

0.07 

T- 9 

1200 

200 

13.7 

1.8 

0.17 

T-10 

1200 

400 

6.8 

3.7 

0.21 

T-ll 

1200 

660 

4.1 

6.1 

0.27 

T-17 

1200 

760 

3.6 

7.0 

0.25 

T-13 

1000 

87 

27.2 

0.92 

0.15 

T-14 

1000 

175 

13.5 

1.85 

0.21 

T-15 

1000 

350 

6.8 

3.7 

0.29 

T-16 

1000 

760 

3.1 

■ 8.0 

0.38 

T-19 

800 

200 

10.0 

2.5 

0.05 

T-20 

800 

400 

5.0 

5.0 

0.08 


Flow Rates: H 2 =200 

(std. cc./rain.) C 0 2 0 2 Q0 



Fig. 52 - Anatase Morphology in Titania Powders, 59,OOOX 


107 



Fig. 53 - Typical Micrograph of Rutile Powder, Run T-8, 40,000X 
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Fig. 54 - Typical. Micrograph of Rutile Pwwder, Run T-6, 40,000X 


109 



110 



Fig. 56 - Typical Micrograph of Rutile Powder, Run T-ll, 40,000X 
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CUMULATIVE FREQUENCY (%) 






Fig. 58 - Effect of Temperature and Residence Time 
on Particle Size of Vapor-Grown Rutile 
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Fig. 59 ~ Volume Growth Rate versus Temperature for Titania Runs 


decrease in particle size with increasing temperature, all other system 
conditions being fixed. The decrease in growth rate below 1000°C indi- 
cated that, in addition to the steadily diminishing nucleation rate, 
reaction and/or condensation kinetics were limiting growth. Few nuclei 
were formed, and their size was very small. 

To test the validity of the assumption that nucleation was 
depleting the gas stream and inhibiting growth, a technique was required 
which would reduce the nucleation rate at a fixed temperature. One 
factor which may affect nucleation is the position of the injector tip. 
Smoke tube studies have been reported 64 which show considerable turbu- 
lence generated as reactant gases exit the injector. This turbulence 
enhances gas mixing and heat transfer and therefore it follows that 
inducing such turbulence near the hot zone may also promote nucleation 
kinetics. Injecting the gases further upstream in a cooler region may 
limit turbulence-induced nucleation by establishing laminar flow into 
and through the hot zone. If nucleation rate can be effectively reduced 
in this manner, a particle size increase would be anticipated i 

Three additional rutile runs were made at 1200°C with the 
injector shortened 8 inches to inject the reactants 10 inches from the 
isothermal zone. Figure 60. Residence times of 0.9, 1.8 and 3.6 seconds 
were employed, corresponding to previous runs T-8, T-9 and T-10. Lineal 
analysis of electron micrographs, Figures 6l to 63 > revealed a consider- 
able size increase over earlier counterparts. Table XXVI. Graphical 
presentation of the two sets of runs. Figure 64, shows the cube of the 
diameter again varies linearly with time. The volume growth rate, G, 
has increased from 2.9 to 19.7 (micron 3 x 10 3 /second) . This nearly 
sevenfold increase in growth rate indicates the importance of injection 
temperature in controlling nucleation and growth. Having confirmed the 
effect of nucleation rate on particle growth, and having negated to some 
extent the nucleation-depletion of the gas stream, a measurably broader 
range of particle sizes was produced. 


Table XXVI - Particle Size Variation with Injector Length 


Run No. 

Residence Time 
(seconds) 

Injection Distance 
from Hot Zone 
(inches) 

Average Particle 
Size 
(micron) 

T- 8 

0.9 

2 

0.07 

T-24 

0.9 

10 

0.26 

T- 9 

1.8 

2 

0.17 

T-21 

1.8 

10 

0.33 

T-10 

3.7 

2 

0.21 

T-25 

3.6 

10 

0.42 
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Fig. 6l - Typical Micrograph of Rutile Powder, Run T-24, 

4o,ooox 
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Fig. 62 - Typical Micrograph of Rutile Powder, Run T-21, 

4o,ooox 
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Fig. 63 - Typical Micrograph of Rutile Powder, Run T-25, 
23,OOOX 
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RESIDENCE TIME (seconds) 



Fig. 64 - Effect of Residence Time and Injector Length 
on Particle Size of Rutile 
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C. ZINC ORTHOT ITANATE (Zn 2 Ti0 4 ) POWDER PRODUCTION 
AND CHARACTERIZATION 

1. Using Zinc Metal as the Zinc Vapor Species 

Zinc orthotitanate powder was produced from homogeneous nuc- 
leation of gases using vaporization of metallic zinc as the zinc vapor 
species. The Zn/Ti ratio shown in Table XXVII is the ratio of the 
number of moles of Zn to number of moles of Ti introduced into the 
system during the run. For stoichiometric Zn 2 Ti0 4 , the Zn/Ti ratio 
should be 2.0. The overall reaction was 

2ZnCl 2 + TiCU + 4H 2 + 4C0 2 * Zn 2 Ti0 4 + 8HC1 + 4C0 . ( 89 ) 


a. Zinc Vaporization - To establish the conditions for a con- 
trolled zinc evaporation rate, the system pressure was maintained at 

50 TmnWg for all r uns - The zinc temperature was varied from 500 to 600°C 
The recorded temperature was approximately 50°C lower than the actual 
temperature . This temperature difference was indicated by the observed 
melting of the metal which occurred at 370 ± 11°C for the runs recorded; 
zinc melts at 420°C . Although the temperature was low, temperature 
measurement gave a relative and reproducible control for establishing 
zinc evaporation conditions. Zinc evaporation rates at 50 mmHg system 
pressure is plotted versus temperature in Fig. 65 . Because of gas flow 
rate variation it was difficult to establish a relation between evapora- 
tion rate and temperature ; although in general, evaporation rate in- 
creased with temperature. 

Several factors may be responsible for poor control of zinc 
evaporation. The first and most probable factor was variation in zinc 
partial pressure. Even though the system pressure was maintained con- 
stant, variation of other gas flow rates caused a variation in zinc 
partial pressure. Increasing gas flow rates would decrease zinc partial 
pressure and increase zinc evaporation rate. This was indicated in runs 
2, 6, and 7 where zinc temperatures were fairly close and zinc evapora- 
tion seemed to be controlled by zinc partial pressure. The second 
factor was reduction of surface area because of powders settling on the 
evaporation surface. This could not be controlled and the magnitude of 
the effect could not be determined. Later runs were made using zinc 
chloride as the zinc vapor species. 

b. Powder Characteristics - In the majority of the rims, 
powder was deposited in three distinct zones in the reaction tube above 
the hot zone. The zones were the top two inches of the reaction tube, 
Zone 1; approximately 2-4 inches below the top of the reaction tube, 

Zone 2; and from approximately 4 inches below the top of the reaction 
zone to the hot zone, Zone 3 (see Fig. 66 for sketch of zones). 

According to visual evidence in most runs Zone 1 contained mainly 
metallic zinc which was powdery and dark gray. Zone 2 contained large 
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Table XXVII - System Conditions and Phase Analysis for Runs Using 
Metallic Zinc as the Zinc Vapor Species 





Gas Flow Rates 



Ron No 

. Reaction 

Zone Clp to Ti 

H, 

C0o 

Zinc Evaporation Zinc Temp. Zn/Ti 

Zn 2 Ti0 4 

Temp. 

(°c) 

(cc/min at 
25 psig) 

(cc/min 

at (cc/min at 

Rate 

(°C) Ratio 

5 psig) 

5 psig) 

(g/min) 

1 

1200 

85 

168 

336 

2.0 

600 7-8 

2 

1200 

100 

200 

400 

0.40 

575 1.1 

3 

1400 

00 

200 

400 

0.16 

500 0.6 

4 

1040 

27 

570 

111 

0.12 

540 1.3 

5 

1100 

100 

500 

1000 


600 

6 

1300 

95 

300 

600 

0.75 

560 2.6 

7 

1350 

100 

300 

1600 

1.90 

570 5-6 

8 

1300 

100 

2000 

1600 

O.67 

520 1.9 


System 


Phase Analysis and Color 

of Powder 



Pressure 

Zone 1 


Zone 2 

Zone 3 


(nmiHg) 

Phases Color 

Phases Color Phases 

Color 

1 

50 

Not Analyzed 


Not Analyzed 

(1) Zn 2 Ti0 4 

(2) Ti0 2 

(3) ZnO 

light gray 


2 

. 50 

Not Analyzed 


Ti0 2 white 

(1) ZnO 

white 





(2) Ti0 2 

(3) Zn 2 Ti0 4 



3 

50 

Not Analyzed 


Not Analyzed 

(1) Ti0 2 

(2) ZnO 

white 






(3) Zn 2 Ti0 4 Trace 

(4) ZnTi0 3 Trace 

4 

50 

Not Analyzed 


Not Analyzed 

(1) Ti0 2 

(2) ZnTi0 3 

light gray 

5 

50 

(1) Zn dark 

gray 

(l) Zn dark 

gray (l) Ti0 2 

light gray 


(2) ZnCl 2 

(3) tio 2 


(2) ZnCl 2 

(2) ZnO Trace 



6 

50 

Not Analyzed 


Not Analyzed 

(1) ZnO 

(2) Zn 2 Ti0 4 

(3) Ti0 2 

(4) Zn Trace 

light gray 






(5) ZnTi0 3 Trace 

7 

50 

Not Analyzed 


(1) ZnO gray 

(1 ) tio 2 

light gray 



(2) Zn?Ti0 4 

(2) ZnO 






(3) 110 2 

(3) ZnTi0 3 . 


8 

50 

Not Analyzed 


Not Analyzed 

(1) TiOa 

light gray 




(2) ZnO 

(3) ZnTiOg 




* All 

zinc evaporated due to variation in system pressure. 
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Zinc Evaporation Rate, g/min 



Fig. 65 - Zinc Evaporation Rate vs Temperature at 50 mmHg 
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Fig. 66 - Schematic of Powder Deposition Zones 


12b 


amounts of metallic zinc which condensed on the walls and flowed down 
the tube to the point that vaporization occurred again. The major 
oxides were found in Zone 3 and were evidenced by a light gray to white 
color. In no run was Zn 2 Ti0 4 the only oxide formed and in all runs Ti0 2 
and ZnO were found mixed with the Zn 2 Ti0 4 . The results of x-ray dif- 
fractometer phase analysis (see Appendix F) as well as material color 
for the different portions of the tube are listed in Table XXVII. The 
indicated phases are listed in order of peak height to indicate relative 
amounts of phases. 

In four of the runs Zn 2 Ti0 4 was formed and in one run (#1) 
it was found to be the most abundant phase in Zone 3* Figures 67 and 
68 are electron micrographs of run Zn 2 Ti0 4 -l, Zone 3* Although the 
characteristic morphology of Zn 2 Ti0 4 is not known and cannot be deter- 
mined from the mixed phases in the micrographs, the micrographs do 
show a particle size in the 0.1 to l.Op. regions. The characteristic 
hexagonal morphology shown by rutile in earlier micrographs is apparent. 

Although other oxide phases were formed besides Zn 2 Ti0 4 , its 
presence indicates that reaction was proceeding, at least to some degree. 
The most probable explanation is lack of mixing of the Zn vapor with 
the gases flawing through the central injection nozzle. The Reynolds 
number calculated for run Zn 2 Ti0 4 -8 for gas flow through the hot zone 
was 60. This indicates laminar flow and suggests a concentration 
gradient across the tube, with a high Zn concentration near the outside 
and a high TiCl 4 concentration in the center. 


2. Using Zinc Chloride as the Zinc Vapor Species 

Formation of Zn 2 Ti0 4 by homogeneous nucleation was achieved 
using ZnCl 2 as the zinc vapor species. The chloride was formed by 
chlorination of ZnO and Zn 2 Ti0 4 . In the former instance, zinc chloride 
was produced by the controlled chlorination of zinc oxide in the internal 
chlorinator below the reaction zone; simultaneously titanium sponge was 
chlorinated at 450°C and transferred at 350°C to the reaction chamber. 

The second method of halide generation involved the controlled chlorin- 
ation of zinc orthotitanate in the internal chlorinator by the following 
reaction. 

Zn 2 Ti0 4 + 4C1 2 + 4C0 * 2ZnCl 2 + TiCl 4 + 4C0 2 (90) 

The feed stock was prepared by solid state reaction between stoichio- 
metric amounts of ZnO and Ti0 2 (anatase). The reactants were fired to 
925°C for 24 hours with only 3 per cent ZnO remaining unreacted. 

System conditions for ZnO source runs are shown in Table XXVIII. 
Flow gauge pressures were as indicated or 0 psig. Conditions for 
orthotitanate sources are shown in Table XXIX. 
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Fig. 68 - Electron Photomicrograph, Run Zn 2 Ti0 4 -l, Zone 3, 82,OOOX 
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Table XXVHI - Syatsa Conditions and Phase Analysis for Runs Using Oxide to Produce Zinc Chloride 




Reaction Zone 





Gas Flow Rates 




Material 

Run So. 

Temperature 

Cc) 

Cl? to ZnO 

a? to Tio? 

00 to ZnO 

CO? 


Ha 

Flow Gage 
Pressure 

(psig) 

Flow 

Rate 

(cc/min) 

Flow Gage 
Pressure 

(psig) 

Flow 

Rate 

(cc/min) 

Flow Gage Flow 

Pressure Rate 

(PSIG) (cc/min) 

Flow Gage 
Pressure 
(PSIG) 

Flow 

Rate 

(cc/min) 

Flow Gage Flow 

Pressure Rate 

(PSIG) (cc/min) 

ZnO 

n 

1300 

25 

40 

25 

0 

5 80 

5 

6 

5 80 

ZnsTio* 

12 

1300 

25 

50 

25 

50 

$ 100 

5 

100 

5 200 

ZnaTiO* 

13 

1300 

25 

50 

25 

50 

5 100 

5 

100 

5 200 

ZnaHO« 

14 

1430 

0 

45 

0 

U5 

0 100 

0 

100 

0 100 

Zn,Tio« 

15 

1400 

0 

38 

0 

38 

0 90 

0 

100 

0 130 

zn*no 4 

16 

1480 

25 

50 

25 

50 

0 17 

0 

17 

0 233 

ZnaT10 4 

17 

1450 

25 

40 

25 

40 

0 100 

0 

100 

0 200 

ZnaTi0 4 

18 

1700 

0 

45 

0 

«»5 

0 100 

0 

100 

0 100 



^rstem 
Pressure 
(mm Hg) 


ZnO Chlorination 



Phase Analysis' 





(*) 



Zone #1 

Zone #2 


Zor.e #3 

ZnO 

U 

210 


46 



Not Analysed 

hot Analyzed 

(1) ZnO 

(2) ZnCl 2 4Za(OH), 

Zi\bT10 4 

12 

100 


53 



ZnCla 

Hot Analyzed 

Nothing Analyzed 

Zn^TlO^ 

13 

200 


61 



ZnCla 

ZnCla 


TlOu 

ZnaTiO* 

14 

100 


70 



ZnCla 

ZnCla 


Tip* 

ZnjTiO* 

15 

90 


101 



ZnCla 

ZnCly 


TiQa 

ZnaT10 4 

l£ 

200 


48 



(1) ZnCla 

(2) TiQ» 

(1) ZnCla 

(2) TlOa 


(1) TlOa 

(2) ZnO 

(3) ZnTiOs 

ZnaT10 4 

17 

300 


57 



(1) ZnCla 

(2) TiQa 

(1) TiOa 

(2) ZnTi0 3 

(3) Trace ZnO 

(4) Trace ZnaTi0 4 
or ZnaTijOo 
or both 

(1) TiCfe 

(2) ZnTiOs 

(3) ZnO 

ZoaTlOe 

IB 

200 


99 



(1) Zn 

(2) ZnCla 

(3) WQ* 

(1) Zn 

(2) ZnCla 

(3) TiQa 


(1) ZnO 

(2) TiQa 

(3) ZnTi0 3 

(4) Trace ZnaTiO 
or ZneTi 3 0a 
or both 



Table XXVIII - (Continued) 


Run No. 

Furnace 

Temp. 

(°C) 

CI 2 to ZnO 
(cc/min) 

CI 2 to Ti 
(cc/min) 

CO C0 2 

(cc/min) (cc/min) 

h 2 

(cc/min ) 

System Pressure 
(mmHg) 

ZO 1 

SC 

1495 

45 

45 .. 

100 

100 

200 

150 

2 

SC 

1450 

100 

50 

100 

50 

200 

350 

3 

SC 

1400 

100 

200 

100 

300 

400 

150 

4 

IF 

1400 

100 

200 

100 

300 

400 

150 

5 

IF 

1300 

200 

200 

- 

200 

4oo 

380 

6 

IF 

1400 

75 

50 

- 

100 

200 

50 

7 

IF 

1250 

100 

200 

200 

600 

600 

200 

8 

IF 

1150 

200 

75 

300 

1200 

600 

150 


Phase Analysis 



Zone 1 

Zone 2 

Zone 3 

ZO 1 

ZnO 

Zn 2 Ti0 4 (tr) 

ZnCl 2 

No Powder 

2 

ZnO, ZnCl 2 , Zn 2 Ti0 4 (tr) 

ZnCl 2 , Zn 

No Powder 

3 

ZnCl 2 , Zn 2 Ti0 4 (tr) 

ZnO, Zn 2 Ti0 4 

No Powder 

4 

ZnO, TiCl 3 

ZnO, ZnCl 4 -Zn(0H) 2 
TiCl 3 

( ZnO, TiCl 3 , Ti0 2 
\ ZnCl 2 -Zn(0H) 2 
< Zn 2 Ti0 4 (tr) 

| Zn 2 Ti 3 0 3 (tr) 

1 ZnTi0 3 (tr) 

5 

Not Analyzed 

ZnO, ZnTi0 4 , TiCl 3 , ZnTi0 3 
ZnCl2-Zn(0H) 2 


6 

ZnO, ZnCl 2 Zn(0H) 2 , TiCl 3 

ZnO 

No Powder 

7 

Not Analyzed 

Not Analyzed 

Ti0 2 

8 

Not Analyzed 

Not Analyzed 

ZnO 
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Table XXIX - System Conditions and Phase Analysis for Runs Using 
Zinc Orthotitanate to Produce Zinc Chloride 


Run 

Temp. 

(°c) 

Cl 2 to 

Zn 2 Ti0 4 CO C0 2 H a 

(cc/min) (cc/min) (cc/min) (cc/min) 

System 

Pressure 

(mmHg) 

0 2 

( cc/min) 

Furnace 

Type 

Zone 

1 

Zone 

2 

Zone 

3 

ZOT 9 

1500 

400 

400 

200 

400 

100 


IF 

NA 

No Powder 

Al 2 Ti0 s 
Ti0 2 ; 

Z-rioTiO 4 

ZnTi0 3 (tr) 

10 

1450 

200 

200 

100 

200 

■ 20 

” 

MF 


ZnCl . 

TiCl 3 

ZnO 

Zn 2 Ti0 4 (tr) 

No 

Powder 

11 

1450 

100 

300 

500 

300 

140 

“ 

MF 


ZnCl 2 

ZnO 

Ti0 2 

ZnTi0 4 (tr) 

Ti0 2 

Zn 2 Ti0 4 

ZnO 

12 

1200 

400 

400 

400 

400 

100 

- 

IF 


No 

Powder 

Ti0 a 

Al 2 Ti0 s 

Zn 2 Ti0 4 

13 

1300 

400 

400 

200 

600 

380 

150 

IF 


Ti0 2 

Al 2 Ti0 s 

TiCl 3 

Ti0 2 

Al 2 Ti0 5 

l4 

1150 

400 

400 

400 

400 

100 

- 

IF 


No 

Powder 

Al 2 Ti0 5 

Zn 2 Ti0 4 

Ti0 2 

15 

1250 

400 

100 

600 

1000 

100 

- 

IF 

ZnCl 2 
Ti0 2 (tr) 

No 

Powder 

Al 2 Ti0 s 

Zn 2 Ti0 4 

16 

1300 

400 

400 

200 

400 

380 

‘ 

IF 

NA 

NA 

Ti0 2 

ZnTi0 3 

Al 2 Ti0 5 

Zn 2 TisP a 

17 

1300 

100 

100 

100 

200 

200 

- 

MF 

NA 

ZnCl 2 

Ti0 2 

Zn 2 Ti0 4 (tr) 




Powder was typically deposited in three zones . Zone 1 was the 
top six inches of the reactor, Zone 2 extended to 12 inches below the 
top of the reactor, and Zone 3 continued to 25 inches from the top. 
Regions listed as "not analyzed" were not barren but contained hydro- 
scopic material identified as zinc chloride. Zinc orthotitanate was 
formed in every run as a minor phase, except in the case of predicted 
mixing conditions which produced 15 weight per cent. 

Figures 69-72 show scanning electron micrographs of sample 

ZOT-12. 

D. CALCIUM TUNGSTATE (CaW0 4 ) POWDER 

Calcium tungstate powder was produced by homogeneous nucleation 
using the chlorination of CaO and W to provide the halide gases . The 
conditions and phase identification for products of formation are listed 
in Table XXX. Figures 73-75 show scanning electron micrographs of the 
products. 


Table XXX - System Conditions and Phase Analysis of CaW0 4 Ruin 


Temperature 

°C 


1500 

Cl 2 to CaO 

cc/min 


100 

Cl 2 to W 

cc/min 


200 

CO to CaO 

cc/min 


100 

C0 2 

cc/min 


700 

h 2 

cc/min 


700 

System Pressure 

mmHg 


100 

Phase Analysis 

Collector 

Zone 1 

Zone 2 


CaO 

CaO 

CaO 


CaW0 4 

CaC0 3 



CaCl 2 

CaW0 4 



Reactor 

CaO 


E. LANTHANUM OXIDE (La 2 0 3 ) POWDER 

Lanthanum oxide powders were produced over the range of l650°C to 
1750°C in the molybdenum element reactor. The furnace heated both the 
reactor and the chlorinator. In runs 1-4, chlorine was passed over 
pure La 2 0 3 ; in subsequent runs carbon was mixed with La 2 0 3 to promote 
reduction and chlorination. With the carbon-free source materials, 
only a trace of La 2 0 3 was detected in a quantity of LaCl 3 ; thus indi- 
cating the need for stronger reduction conditions by the addition of CO. 
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Fig. 71 - Scanning Electron Micrograph, Sample ZOT-12 

5000X 



Fig. 72 - Scanning Electron Micrograph, Sample ZOT-12 

10,000X 



Fig. 73 - Scanning Electron Micrograph, Sample CW-1 

1000X 


Fig. 7^ - Scanning Electron Micrograph, Sample CW-1 

3000X 
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Fig. 75 - Scanning Electron Micrograph, Sample CW-1 

5000X 
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Table XXXI - System Conditions and Phase Analysis for 1*303 Runs 


Run 

Temp. 

°C 

Cl2 

cc/min 

CO 

CDs 

cc/min 

H- 

cc/min 

Carbon 

System 

Pressure 

mrallg 

xs 

Clo 

Zone 

1 

Zone 

2 

1 

1750 

45 


100 

100 

No 

150 

• 


LaOCl 
LaCl 3 
LagOs (tr) 

2 

1750 

100 

- 

500 

500 

Ho 

300 

- 


LaCls *7H 20 
LaOOH 

3 

1700 

100 


300 

300 

No 

50 



LaOCl 
LaCl 3 ‘7H20 
LaOOH 
1*^03 (tr) 

4 

1700 

100 

- 

800 

600 

No 

500 

- 

LaOCl 

LaAlOs 

LaOCl 

LaAlOs 

5 

1700 

300 

100 

- 

300 

No 

30 

- 

LaOCl 

LaCla’THgO 

LaOCl 

LaCl 3 -7H20 

6 

1700 

300 

400 

" 

300 

No 

250 

“ 

LaOCl 

L&CI 3 

La^a 

LaOCl 

LaCl 3 

La^Ps 

7 

1750 

300 

• 

500 

300 

Yes 

50 

• 

L&CI 3 

LaOCl 

LagOa 

LaCl 

LaOCl 

LagOs 

8 

1750 

300 


500 

300 

Yes 

100 


LaOCl 

LaCls 

LaA10 3 

LaOCl 

LaCl 3 

LaAlOs 

I*20 3 

9 

1650 

300 

300 

500 

500 

Yes 

50 

- 

LaOCl 

LaOCl 

10 

1700 

300 

200 


300 

Yes 

50 


LaOCl 
LaaOs (tr) 

LaOCl 
LaOOH 
LaA10 3 
LagOs (tr) 

11 

1750 

100 

100 

700 

700 

No 

100 

• 

LaOCl 

La^Os 

LaCls* 7 H 2 O 

LaOCl 

LaOCH 

12 

1750 

200 

100 

300 

300 

Yes 

100 


LaOCl 

LaOOH 

LaAlOs 

LaCl 3 *7HaO 

LaOCl 

LaOOH 

LaAlOs 

LaCl 3 *7H20 

13 

1750 

150 

300 

300 

300 

No 

150 


LaOCl 

LeCls-7HaO 

LaAlO 

LaOOH 

LaOCl 
LaCls -7H20 
LaAlOs 

14 

1750 

50 

300 

300 

300 

No 

150 

' 

NA 

LaOCl 
LaCl 3 -7H20 
0*203 (tr) 
LaA10 3 

15 

1750 

150 

500 

300 

300 

No 

380 

50 

LaOCl 

I*2p3 

LaCl 3 *7HgO (tr) 

LaOCl 

1*20 3 

LaCl 3 .7H 2 0(tr) 

16 

1700 

100 

300 

300 

300 

No 

380 

50 

LaOCl 

LaOCl 

17 

1750 

200 

100 

* 

300 

No 

100 

50 

LaCls-TOsO 

LaOCl 

LaOOH 

I*Cl 3 *7H 2 0(tr) 

LaOCl' 

LaOOH - 

18 

1700 

200 

200 

" 

500 

No 

380 

50 

LaCls -71120 
LaOCl 

1*203 (tr) 

LaOCl 

19 

1750 

150 

400 

200 

500 

Ho 

250 

- 

LaOCl 

LaOCl 

20 

1750 

220 

* 

500 

300 

Yes 

480 

“ 

LoC1 3 -7H20 

LaOCl 

LaCls -7H20 
La^Os (tr) 

21 

1700 

220 

200 

200 

300 

No 

100 

- 

UC1 3 

LaOOH 

LaCl 3 

LaOOH 

22 

1750 

100 

- 

300 

300 

Yeo 

380 

- 

LaOCl 

(tr) 

LaOCl 

La .703 
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Fig. 76 - Scanning Electron Micrograph, Sample LA.-16 

500X 



Fig. 77 - Scanning Electron Micrograph, Sample LA-l6 

3000X 





Fig. 78 - Scanning Electron Micrograph, Sample LA-20 

500X 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The following conclusions can be made concerning this effort in the 
development and preparation of space-stable pigments for thermal control 
coatings: 

1. Aluminum Oxide Powder 

a. For constant reactant gas inputs and reaction zone 
temperature, increasing reaction zone pressure 
increases the average particle size of alumina 
powder . 

b. Increasing reactant gas inputs and increasing 
system pressure to maintain constant gas veloci- 
ties at constant temperature increases the aver- 
age particle size of alumina powder. 

c . The particle size range for the various alumina 
runs was 0.1 to 0 . 50 p. 


2. Titanium Dioxide 

a. Systematic variations of particle size were con- 
trolled with hot zone temper attires, residence 
time, and injector design. Sizes from 0.06 - 0.42 
micron were produced, with distributions as narrow 
as any yet reported. 

b. A constant volume growth rate of particles at a 
specific temperature was observed. Variation of 
growth rate with temperature reflected the effect 
of changing super -saturation and nucleation rate. 
Additional flexibility in controlling nucleation 
was achieved by varying injector length. 


3. Zinc Orthotitanate 

a. Improvements in the yield of zinc orthotitanate 
have been attained by reactor design improvements 
that incorporate maximized jet action. 

b. Gas mixing must be increased to eliminate the 
mixed species. 
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4. Calcium Tungstate 

Calcium tungstate was produced as the major phase. 

5. Lanthanum Oxide 

a. Lanthanum oxide has been produced as a minor 
phase only. 

b . Lanthanum oxychloride kinetics prevents the 
formation of lanthanum oxide in large amounts . 


B . RECOMMENDAT IONS 

Gas mixing remains the major limitation to high yields of multiple 
oxides . A plasma cell should be incorporated into the present reactor 
as an upstream turbulence and heating source. Further reaction system 
improvements should permit pilot production studies at an early date. 
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FREE ENERGY VALUES USED IN THERMODYNAMIC CALCULATIONS 
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Table XXXII - Free Energy Values Used in Thermodynamic Calculations 



T 

•c 

T 

•K 

Aia 3 

AljO, 

Oa 

00 

00a 

HCL 

Ha 

TiCl* 

TiQa 

In 

ZnClz 

ZnO 

ZnaTiO* 


Ceirc>4 

LaOCl (est) 

LaCl 3 

CeCl 2 

VCls 


27 

300 

-136.0a 

-378.00 

0 

-32.80 

-94.25 

-22.75 

0 

-175.90 

-212.00 

0 

-88.45 

-76.10 

-369.60 

-402.60 

-331.80 

-223.63 

-238.30 

-179.65 

-60.00 


127 

400 

-134.85 

-37O.3O 

0 

-35.00 

-94.30 

-22.95 

0 

-170.55 

-207.90 

0 

-84.80 

-74.65 

-362.20 

-393-10 

-323.00 

-218.70 

-233.00 

-175.95 

-52.50 


227 

500 

-133.61 

-362.80 

0 

-37.10 

-94.40 

-23.90 

0 

-167.45 

-203.60 

0 

-81.10 

-71.30 

-351.20 

-384.50 

-314.60 

-213.97 

-227.40 

-172.50 

-45.00 


327 

600 

-132.36 

-355.30 

0 

-39.35 

-94.45 

-23.40 

0 

-164.50 

-199.30 

0 

-78.80 

-68.95 

-342.20 

-376.60 

-306.15 

-209.77 

-222.70 

-169.05 

-40.40 


427 

700 

-131.09 

-347.80 

0 

-41.55 

-94.45 

-23.60 

0 

-161.60 

-194.95 

0 . 

-76.60 

-66.60 

-333.15 

-368.10 

-397.80 

-205-13 

-217.30 

-165.75 

-37.80 


527 

800 

-129.81 

-340.40 

0 

-43.70 

-94.50 

-23.75 

0 

-158.70 

-190.75 

0 

-74.40 

-64.00 

-323-75 

-360.00 

-289.50 

-200.43 

-211.30 

-162.35 

-35.20 


527 

900 

-128.51 

-332.90 

0 

-45.85 

-94.55 

-23.90 

0 

-155.85 

-106.55 

0 

-72.20 

-61. 50 

-31“ -50 

-351-30 

-281.30 

-196.10 

-207.00 

-159.15 

-32.60 


727 

1000 

-127.01 

-325.20 

0 

-47.95 

-94.60 

-24.00 

0 

-152.95 

-182.35 

0 

-70.00 

-59.00 

-305.35 

-342.40 

-273.10 

-191.10 

-200.90 

-155.70 

-30.00 


OH 

1100 

-125.41 

-317.20 

0 

-50.10 

-94. 60 

-24.15 

0 

-150.05 

-178.00 

0 

-66.72 

-56.10 

-295.20 

-334.90 

-265.10 

-186.90 

-195.86 

-152-80 

-27.30 

H 

■pr 

927 

1200 

-123.00 

-309.40 

0 

-52.15 

-94.65 

-24.40 

0 

-147.10 

-173.89 

0 

-63.44 

-53-40 

-285.69 

-326-90 

-256.95 

-182.83 

-191.60 

-150.30 

-24.70 


1027 

1300 

-122.10 

-301.50 

0 

-54.35 

-94.70 

-24.40 

0 

-144.15 

-169.55 

0 

-60.16 

-50.70 

-275.95 

-319-00 

-248.80 

-179-00 

-188.00 

-147.70 

-22.10 


1127 

1400 

-120.57 

-293.30 

0 

-56.25 

-94.75 

-24.65 

0 

-141.15 

-165.45 

0 

-56.88 

-47.90 

-266.25 

-310.50 

-240.60 

-175.00 

-184.50 

-145.10 

-19.50 


1227 

1500 

-198.95 

-286.10 

0 

-58.40 

-94.75 

-24.70 

0 

-138.15 

-161.20 

0 

-53-60 

-45.10 

-256.40 

-302-50 

-232-35 

-171.13 

-180.90 

-142.40 

-17.00 


1327 

1600 

-117.31 

-278.10 

0 

-60.60 

-94.80 

-24.90 

0 

-135.20 

-157.25 

0 

-51.50 

-42.40 

-246.65 

-294.50 

-224.50 

-167.23 

-177.20 

-139.60 

-14.20 

«:d 

1427 

1700 . 

-115.67 

-270.60 

0 

-62.65 

-94.85 

-25.00 

0 

-132.25 

-153.10 

0 

-49.50 

-29.40 

-236.90 

-286.00 

-216.45 

-163-23 

-173.70 

-137.00 

-U.5O 

O 

w*r« 

1527 

1000 

-111*. 03 

-263.10 

0 

-64.65 

-94.90 

-25.25 

0 

-129.25 

-149.35 

0 

-47.50 

-36. 60 

-227.55 

-278.00 

-209.05 

-159.40 

-170.20 

-134.30 

-9.00 

'O 

try 

a—i 

0 

1627 

1900 

-112.30 

-255.50 

0 

-66.75 

-94.90 

-25.45 

0 

-126.25 

-145.25 

0 

-45.50 

-33.80 

-217.85 

-270.00 

-202.00 

-155.57 

-166.70 

-131-70 

-6.40 

1727 

2000 

-110.73- 

-247.90 

0 

-68.75 

-95.00 

-25.50 

0 

-123.25 

-141.15 

0 

-43.50 

-31.00 

-208.15 

-262.00 

-195.05 

-151.67 

-163.00 

-129.00 

-3.80 
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APPENDIX B 


COMPUTER PROGRAM SOURCE LIST FOR CALCULATION OF 
EQUILIBRIUM CONSTANTS AND FREE ENERGY VALUES 



PAGE BLANK NOT 1 



t- 


FCRTRAN IV G LEVEL It MOD 3 MAIN DATE = 69071 11/29/22 


0001 DIMENSION I DLGR ( 20 ,20 ) , I DLGP ( 20 , 20 ) , I R ( 1 5) , I P ( 15 ) 

0002 DIMENSION I A ( 1 5 ) , I B ( l 5 > , IC ( 1 5 ) , MA ( 1 5 ) , MB ( 1 5 ) , MC ( 1 5 ) 

000 3 DIMENSION JA( 1'5) , JC( 15 ) ,KAtl5) ,KB(15) ,KC( 151 

0004 999 FORMAT ( 15,1515) . _ 

0005 998 FORMAT < 3 X » 3A4, 91 7) 

0006 997 FORMAT! • DELTA G IS* » 1 8 , • CALORIES AT*, 13, ’00 DEGREES KELVIN AND T 

I HE EQUILIBRIUM CONSTANT IS',F12.4) 

0007 995 -FORMAT (.*0 THE FOLLOWING DATA IS FOR THE PF ACT ION *) 

0008 994 F0RMAT! , OPEACTANTS , »4X,I?,3A4,I2,3A4, I 2 » 3 A 4*, 12 ',3 44, I 2 , 3 A 4 ,12, 3 A 4 ) 

0009 993 FORMAT! 'OPRGDUCTS' ,5X, 12 ,3A4, 12, 344, I ?,3A4, 12 ,3A4 , I 2, 3A4, I? ,3A4// ) 

001 0 992 FORMAT! *1 THE DATA USED IN THFSE C ALCUi.AT f CNS IS ') 

0011 991 FORMAT! *0* ,2X,3A4, 917) 

0012 READ! 5,999) IRCTN,! IR ( I ), 1 = 1, IRCTN) 

0013 RE AD (5, 999) I PRDT , { I P ( I ) ,1=1, IPRDT ) 

0014 R E AD ! 5 , 998 ) < ( I A ( I ) , I B ( I ) , I C ( 1 !TT I DLGR ( I » M ) , M= 3~, 1 HYMA ! If, MB ( I ) , MC! 

II),! I DLGR ( I ,M),M= 12,20 )) ,1=1 , IRCTN) 

0015 READ"! 5,598) ( ( J A ( I ) , J B ( I ) » JC ( I ) , ( I DLGP! I , m ) , M= 3 , 1 1 ) , K A ( I ) , KP ! I ) , KC ( 

II ), t IDLGP! I »M) »M=12,20) ) ,1=1 , IPRDT) 

"'001 6 90 WRITE(6,992) ' " 

0017 92 WRITF!6 ,091) < { 14 < I ), IB! I ) ,IC !T) , ( IPLGR ( I ,M ) ,M= 3, 1 1 > ,MA! I ) ,MR( I },MC 

1! I), I IDLGKU ,M),M=12,20) ) ,1 = 1, IRCTN) 

0018 91 WRI TE (6,991 ) l ( JA ( I ) ,JB(I ) , JC ( I) , ! IDLGP (I ,M) ,M = 3, II ) ,KA( I ) ,KB! I ) , KC 

1 ( I ) , ( IDLGP! I ,M) ,M=I2,20) ) , 1 = 1 ,'IPROT > 

0019 95 WRITE(6,995) 

0C20 94 WR I TE (6,994 ) ( ( I R ( I ) , I A! I ) , I B ( I) , IC ! I ) ) , 1= I , IRCTN) 

002 1 93 WRI TE(6 ,993 ) ( ( IP ( I I , JA! I ) , JB ( I ) , JC. ( I ) ) , 1= 1 ♦ I PPCT ) 

0022 99 P096M=3 ,20 

002 3 A =0 

0024 D098K=1 , IRCTN 

0025 98 A=A+IR(K )*IDLGR(K,M) 


0026 B =0 

0027 D097K= 1 , I PRDT 


0028 

0029 

97 B=B+IP!K)«IDLGP(K,M) 
I DLG = B-A 


- 

0030 

0031 

DLG= IDLG 
R=1 .987 



0032 

DFNM=R*M*100 



0033 

X=-DLG/DENM 



0034 

C=FXP (X ) 



00^5 

96 WRITE(6, 997) IDLG, M,C 



00 3 6 

STOP 



00 3 7 

C T 
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Preceding page blank 

// (3000,0), CLASSIC 
//ST1 FyFC PROC = FORTR(IM, 

// RFG!on.go=?5?k,timf.go=(o,50) , 

// P A R M . G 0 = • SI 7 F = 1 5 0 K ' 

//cmp. sys in nn * 

DIMENSION C(70,?0), CTFMP (70,70), RCTFMP ( 70, 70 ) , INDEXRI70), 

1L TEMPI 70) , LPR ( 70) , I. PC (70) ,RNV( 70, 70 ) , 0 ( 70 ),RLNK ( 70) ,RK( 70) , 

7GAMMA ( 70 ) , IVFCTR (70 ) ,PAPP t 70 > , IFLFM ( 70 ) 

COMMON NC0LM1 , RNV , C , DFL T AZ , INOFXR,NCOL 

FPSLOMr-1 .OF-OB 

1 RF AD 1 000, MROH, MCOL , UMAX, MCASF _ 

1000 FORMAT (4110) 

IF ( FOF , 60 ) 10,11 

10 CALL FXIT 

11 NC0LM1=MC.ol-1 

RF AD 3000, ( IFLFM(I) , I =1 ,MCOLMl ) 

3000 FORMAT (PA 10) 

RF AD )001 , ( ( C< 1 , J ) , J=1 ,MROU ) , J=1 ,NC0LM1 ) 

1001 FORMAT ( PF10. 1 ) 

RF AD 1001 , (GAMMA! I) , 1=1 ,MROW ) 

DO 500 ICASF=1 , MCASF 
RFAD 1007, P,T 
100? FORMAT ( 7 F 1 0 . 1 ) 

RFAD TOO), (G( I) ,1=1, MROH) " ~ ” 

RFAD 1001 , (C ( I ,NCOL ) , 1 =1 ,NROW ) 

PRINT 1099, P,T 

1999 FORMAT ( 1H1 , RHPRFSSORE , E20. 8, 10X , 11 HTFMPFRATIJRF ,F?0.8) 

PRINT 7000 

2000 FORMAT ( 10X, 1 PHCOMPOS IT ION MATRIX) 

PRINT 3001 ,( !FLFM( I ) , 1=) ,NC0LM1 ) 

3001 FORMAT ( P A 1 0 ) 

DO 1 7 1=1 , MROW 

12 PRINT 2001, (C ( I , J ) , J=1,NC0LM1) 

2001 FORMA 1 (IX, 1 OF 10. 7) 

PRINT 2002 

7002 FORMAT I///1X.70H GIBBS FRFF ENERGIES ) 

PRINT 2003, (G(I ), 1 = 1, MROW) 

700? FORMAT (1X,BF14.6) 

PRINT 2004 

7004 FORMA! <///lX,17H PHASF PARAMETERS ) 

PRINT 7001 , (GAMMA! I ) , I =1 ,NROW) 

PRINT 2005 

2005 FORMAT (//1X.30H FIRST ESTIMATE OF COMPOSITION ) 

PRINT 7003, (C! I ,NCOL ) , 1=1 ,NROW ) 

PRINT 7 00 A 

2006 FORMAT! // 1 3X , 7HSPFC I FS , 9X , 1 1 HCOM POS I T lOM , 5X , 70HF0II I L I RR IIJM CONST AN 
IT ,5X, 1 6HPART IAL PR F S SDR F , 10X , 8HC VFCTOR) 

DO INITIALIZATION 


I BAS I S = 1 
1C0NV =0 
I TF R = 0 

70 IF(ICONV) 71,21,2? " 

22 TOTP=0 

DO IP I = I, MROW 

PARPI I )=A*C( I , N C 0 L ) *GA M M A ( I ) 

IP TOTP=TOTP + PARP ( I ) 

PRINT 70??, ( I,C ( I ,MCOL ) ,RK( I ), PARPI I), IVFCTRI I ) , 1=1, MROW) 

702? FORMA! ( 1 X, II 9, 7F70.P.F74.P, I?0) 

PRINT 7773 , A , TOT P 

???? FORMA! ( / 1 OX , 1 7HVAL0F OF A = , F 70. 7 , 1 OX , 1 4HT0TAL PRESSURE ,F?0.7) 
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PRINT 2222, ITER 

2222 FORMA 1 (/ 10X, 33HTHESF RESULTS WFRE OBTAINFD AFTFR IA.12H 1TERATIO 
INS ) 

CO TO 500 
21 ICONV =1 

IF I ITFR-JTMAX) 25,23,2? " "" — 

23 PRINT 2023, ITMAX 

2023 FORMAT (1 OX, 15HITFR AT ION LIMIT, 15, 10H RFACHFO ) 

GO TO 22 
25 J TER = ITER+1 

. IF ( IliASIS ) 199, 199, 49 

C STORE THE COMPOSITION MATRIX AND THE VFCTOR FSTIMATFS OF COMPOSITION 

C IN THE ARRAY CTFMP 

C 

49 00 40 J=1 , MROW 
00 41 1=1, NCOL 

41 CTFMPI J, I )=C< J, I ) 

40 INDEXR ( J ) = J 
C 

C CONSTRUCTION OF OPTIMUM BASIS 

C 

NROW 1 = NROW 

50 IF ( I-NR0W1 ) 69,69,48 
68 PRINT 2068 

2068 FORMAT (1X.37HNO NONSINGULAR B MATRIX CAN BF FOUND ) 

GO TO 500 
69 MAX=:. 

no si k = i , nrow i — — 

IF (CTEMPIMAX.NCOL I-CTFMP (K ,NCOL 1)52,51,51 
52 MAX=K 

51 CONTINUE 

00 56 N= 1 , NCOL M ] 

56 RC TEMP I I ,N)=CTEMP(MAX,N) 

C CHECK FOR LINEAR OEPFNOENCF 

IF ( I - 1 ) 54,661,54 

54 1 M ) = I -] . 

00 57 N = 1 , I Ml 

RTFRM=0.0 1 

On 58 I SUM= 1 , NC0LM1 

58 RTFRM = RTERM + RCTFMP ( 1 , I5UM )*RCTFMP(M, I SUM ) 

00 59 J = 1 , NCOL Ml 

59 R CTFMPI I,J)=RCTFMp( j ,J)-RTFRM*RC TEMPI M,j) 

_ 57 CONTINUE 

661 RT FRM=0 

DO 60 J = 1 , NCOL M 1 

60 R1ERM=RTFRM+R CTFMPI I , J ) **2 
RTERM=SORTIRTFRM) 

IF IRTERM-EPSLON ) 62,62,61 
C 

C RTERM-E PSLON LFSS THAN OR EQUAL 7FR0 IMPLIFS DFPFNDENCE 
C RTERM-EPSLON P0S1TIVF IMPLIES INDEPENDENCE 
C 

62 DO 63 J=1,NC0L 

T FMP = CTFMP(NROWI , j ) 

CTEMPI NR0W1 ,J)=CTEMP(MAX, J I 

63 CTFMP(MAX,J)=TFMP 

1 1 TEMP= INOFXR (MAX ) 

JNOFXRI MAX ) = INOEXR I NR0U1 ) 
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INDEXR (MRDW1 ) = T I TEMP 
NPOW I = MROW 1 - ] 

Gn Tn so 

6i rn 66? j=i,ngolmi 

662 Rf. TEMPI I , J)=RC TEMPI I , J I/RTERM 
IF (I -M/NX) 66,66,66 
66 DD 66 N=l,NCf>L 
T E P P - C T E M P ( I , N ) 

CTEMPI I , N ) =CTEMP ( MAX ,N ) 

66 C T E M P ( M A X » N ) = T E M P 
. 1 I TEMP = IMOFXR (MAX ) 

INDEXR I MAX ) = INOFXR (I) 

INDEXR I T ) - I ITEMp 
65 ) F (I-NCOLM1) 67,70,67 

67 1 = 1 + 1 
GO TO 50 

. 70 CONTINUE 

CONSTRUCTION OF RFACTIOM COEFFICIENT MATRIX AMO EQUILIBRIUM 
CONSTANT VECTOR 

COMPUTE B INVERSF AMO STORE THE RFSIJLT IN CTEMP 

call I N V 1 (CTEMP, NCOLM 1 ,70,0,LTEMP,IERR,DET,NPIV,PIV,LPR,LPC) 
1 F ( IFRR ) 91,92,91 

91 PRINT 2091, IFRR 

2091 FORMA! ( 1 X , 6H I E RP, , I 5 , 5X , 36HTROUBL F CONSULT DEFINITION OF I ERR 
GO TO 500 

92 DO PO 1 = 1, NCOLM! 

DO HO J= 1 , MROW 

T F M P = 0 ' “ 

DO 81 N= 1 , NCOL M 1 
PI TFMP=lEMP-tC( J,N)*CTFMP(N,I I 
HO RNV( J, I ) =TFMP 
DO 86 1 = 1, NROW 
TFMP=0 

DO 86 ,1 = 1 , NCOLM! 

IN=INDEXR(J) 

86 TFMP=1EMP+RNV( ! , J)*G(]N) - 

R=l. 98726 

RLNKt )) =-1.0/ (R*T)*< GII) -TEMP) 

85 RK ( I > = EXP ( RLNK ( ID 


EQUILIBRIUM PROCEDURE 

199 T E M P = 0 

DO 200 I = ] , N R 0 W 

.200 TEMP = TEMP + GAMMA ( I )*C( I ,NCOL ) 

A = P / 1 E M P 

DO 19 J = 1 , N ROW 
19 1VEC1RI J )=0 

DO 301 J=I , NC0LM1 

IF ( I-INDFXR(J) 1301 ,600,301 

301 CONTINUE 

1 E M P = 0 ‘ 

DO 302 J=1,NC0LM1 
1N=1NDEXR(.D 

IF ( GAMMA ( I N )*RNV( I , J ) ) 307,308,307 
308 T F R M = 0 

GO TO 30? 

307 1 FRM = GAMMA ( IN ) *RNV ( I , J ) =» AL OG ( A$C ( IN ,NCOL ) ) 

302 TFMP= I EMP+TFRM 

IF ( GAMMA (II) 306,309,306 
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C 

C 

C 


309 TFRM=0 

GO TO 255 

306 TFRM=GAMMA( I ) *ALOG ( A*C ( I ,NCOL ) ) 

255 RLNP=-TFMP+TERM 
RO=EXP!RLNQ) 

IFfGAMMA! I ) 1303,320,30? 

303 on 304 J=1,NC0LM1 

IF (RNV(I.J)) 305,304,^05 
305 IM=IMDFXR(J) 

304 .cn F N nm!r CnL, ‘ 0 * 01 ^ C,1M ’ NCPL,,30 4* 30 4.32o 

TREATMENT OF MINOR SPECIFS 
NMULT= ITER/-, 

IF ( ITER-NMULT*4 ) 350,310,350 

FlAp A ri 1 ’ /R0, -1*°' : -04) .300,300,312 

ET AP = C( I , NC.ni ) *RK ( I ) /RO 

DFL TAZ = FTAP-C ( I , NCOL ) 

CALL FUDGE! I) 

C( I *NCOL )=C( I, NCOL J+OFLTAZ 
1 VEC T R ( I ) = ] 

..DO 313 J = 1 , NCOL Ml 

I N= I NOE XR ( J ) ‘ - — 

C( I N ,NCOL ) =C ( IN, NCOL )-RMV( I , J )*DFLTA7 
ICONV=0 4 

GO TO 300 


310 

312 


313 

350 


320 

321 

322 


324 
327 

333 

340 

325 


C 

C 

C 


TREATMENT OF MAJOR SPECIES 

i 1 )^* I,/R0, ~ 1 *° F ' 05> 300 '300,321 

IF (GAMMA! I)) 322,323,3?? 

T EMP=0 

DO 327 J = 1 , NCOL Ml 

IN= I MDEXR ( J) - - - 

cnN?lLfjE P+0AMMA< IN,!< ‘ RNV< in,mcoli 

TEMP=TEMP+GAMMA( I )/C( I , NCOL ) 

OFLTAZ= (RLNk! I ) -RLNP ) /TFMP 
CALL FUDGE! I) 

C ( 1 , NCOL ) =C( I »NCOL ) +DELTAZ 
00 325 J = 1 , NCOL M 1 
JN=1 MDEXR (J) 

C ' IN ’ NC " L ,=C< IN *NCOL l-RNV! I, J)*DFLTAZ 
GO TO 350 


323 

330 

337 

'331 


332 

334 

386 


CONDENSED PHASE PROCEDURE 

IF ( C ( I, NCOL 1-1 . 0E-06 ) 330,330,331 
IF (RLNK(I)-RLNO) 337,331,331 
IVECTRt I ) =4 ’ ' 

GO TO 300 

T FMP = 0 — 

I VECTR ! I ) =3 

DO 334 j = i ,NC0LM1 ' 

IN= INDEXR ( J ) 

IF ( C ( IN, NCOL ) ) 332, ’34, 332 

CrNT = mjE P * r, * MMM ,N,3fiNV< I w. Menu) 

IF ( C ( I , NCOL ) ) 3R6 » 385 , 386 
1FMP=7EMP+ GAMMA! I )/C( I ,N C DL ) 
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INDEXR (NR0W1 )=I I TEMP 
NROW ] =NROWl-l 
GO TD 50 

61 DO 66? J=],NC0LM1 
.. 662 RG.TEMP (I , J ) =RCTFMP ( I , J l/RTPRM 
IF ( I-MAX ) 66,66,66 

66 DD 64 N=1,NCDL 
1 EMP=CTEMP ( I ,N) 

CTFMP( I*N(=CTFMP(MAX,N) 

64 CT EM P (MAX ,N ) =1 FMP 
... ITTFMP= INOFXR (MAX ) 

1 NDE X R ( M A X ) = I NOF XR ( I ) 

JNDE XR ( I )= I I TEMP 

65 IF (J-NCOLMi) 67,70,67 

67 1=1+1 
GO TO 50 

70 CONTINUE 


CwISTect™ RFACTlnN CnFFFrCIFNT MATRIX AND EQUILIBRIUM 
COMPUTE R INVERSE ANO STORE THE P.FSULT IN CTEMP 


] E Rr') ^ 9 9?^ 9 1 * 0 * *- ’ I ^ RR * OFT , NP IV,PIV,LPR,LPC) 

91 PRINT 2091, IFRR 

?09! FORMAT 5 UX,4HIFRR,I5,5X,36HTRO()BI.F CONSULT DEFINITION OF I ERR 

9 ? no 80 1 = 1 , NC.OL M 1 - 

- on 80 J= 1 1 NROW 

DO 81 N=1 ,NC0LM1 

81 T FMP= 1 EMP + C ( J , N ) *c TFMP ( N , I ) 

80 RN V ( J , I ) = TFMP 

00 8 5 1=1 , NROW - 

T E M P = 0 

no 86 j = ] ,MC0LM1 ' " — - : - - 

IN= I NOEXR ( J ) 

86 1 FMP = 1 EMP + RNVI I , J )*G( IN I " "" ~ ■ - 

R=l. 98726 


RLNK( Il=-1 .0/(R*T)*(G(I)-TEMP) 

85 RK ( I ) =EXP ( R L M K (II) 

G EQUILIBRIUM PROCFOURE 

c 

199 TEMP=0 

DO 200 1=] ,NROW 

200 TEMP = TEMP + GAMMA ( I ) =>C ( I , NCOL ) 

A = P/T FM P 

00 19 J= 1 , N ROW 
19 IVECTR( J 1=0 

no 301 J= 1 , NCOL Ml 
IF ( I -IMOFXRf J ) ) 301 ,400 ,301 
301 CONTINUE 

1 F M P = () - - 

no 30? J= 1 , NCOL Ml 

I N= I NOE XR ( «l I 

IF (GAMMA( INI*RNV( I , J » ) 307,308,307 
308 T FRM = 0 

GO TO 30? 

30? I IN, Menu , , 

IF (GAMMA! I)) 306,309,306 
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385 CONTINUE 

J F (1FMP) ?3?,?35,?3? 

335 KM IN = C< I, NCOL ) 

PP 380 J=),NCOLMl 
I N = I ND E XR 1 .1 ) 

I F ( RNV I I , J ) ) 381,380,381 

381 IF (RMIN-Ct IN.NCPL )/ABS(RNV( 1, J))» 380,380,351 
351 RM 1 N = C I I N , NC OL 1 / A B S ( R N V ( I , J ) ) 

380 CONTINUE 

IF (RLMK(I)) 353,300,354 

353 DEL T AZ = -RM I N . 

r,n in 340 

354 OELTAZ=RMIN 
GO 10 340 

300 CONTINUE 

CHECK IF UPDATING OF BASIS IS NECFSSARY 

DO 400 J= 1 , NC0LM1 
I N= 1 NDE XR ( J ) 

IF (RNV( I ,J ) 1410,400,410 
410 IF ( C ( I , NCOL ) -C ( I N , NCOL ) ) 400,400,430 
420 1 R A S I S = 1 
400 CONTINUE 
GO TO 20 
500 CONTINUE 
GO TO 1 
END 

.SUBROUTINE FUOGF(I) 

COMMON NCniMl , RNV,C,nELT#Z , I NDFXR ,NCOL 
DIMENSION INDEXR(TO) , C f 70, 20 ) , RNV ( 70, 20 ) 

IF (C( I, NCOL 1+DFLTA7 1 1,1,2 

1 DEL T AZ = -C< I,NCni.)/2.0 

2 CONTINUE 

DO 3 J = 1 , NC0LM1 

1 N= 1 NDE XR ( J ) 

I F ( C ( IN, NCOL )-RNV( 1 , J 140ELTAZ 1 4,4,3 

4 IF (CON, NCOL)) 5,6,5 
6 DELTA Z =0.0 

GO 10 3 

5 DELTAZ=C<IN,NCOL )/(2.0*RNV< I, J) ) 

3 CONTINUE 
RETURN 
FND 

N-THF ORDER OF A 

SUBROUTINE 1NV1 I A , N , NN , F PS , L TEMP, I ERR , OF T , NP I V, P I V, LPR , L PC ) 
DIMENSION A ( NN , NN ) 

DIMENSION L TEMPI 1 ) ,LPR( 1 ) ,LPC< 1 ) 

C 

C DECK 7046A 
C 

C SUBROUTINES CALLED - NONE 

C THIS ROUTINE INVERTS MATRIX A IN ITS OWN SPACF. IT ALSO COMPUTES THE 
C THF DETERMINANT OF A. 

C 

C THE METHOD IS THF USUAL GAUSSIAN EXCHANGE PROG.FSS. BOTH ROWS AND 
C COLUMNS ARE SEARCHED FOR MAXIMAL PIVOTS. THFRF IS ND UNNECESSARY 
C INTERCHANGING OF ROWS OR COLUMNS, ALL SUCH INTERCHANGES BEING CARRIED 
C OUT AFTER THF EXCHANGE PROCFSS IS COMPLETF. CHAPTFR 1 OF E.L. STIFFLF, 
C INTROPUCT ION 10 NIJMFRICAL MA THFM A T JC S , AC A DFM I C PRF SS, N.Y ., 1 963, SHOULD 
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7046A023 

7046A024 


7046A004 
7046A005 
7046A006 
7046A00 7 
7046A008 
7046AD09 
7046A0 ] 0 
7046A01 1 
7046A01? 
7046A01 3 
7046A014 
7046A01 8 
7046A0 1 6 
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C RE HELPFUL IN FOLLOWING THE CODE. 

C 

C THE CALLING PROGRAM MUST SET A , N , NN , F PS , AMD LTFMP TO- 

C 

c a- the matrix to re inverted 

C NN-THE: NUMRFR DF WORDS OF STORAGE PROVIDED FOR EACH COLUMN 

c OF ARRAY A RY THE CALLING PROGRAM 

C 

C FPS-A NON-NEGATIVE MUMPER WHICH EACH PIVOT IS RFOUIRFD TO EXCEED 

_P ._ IN ABSOLUTE VALUE ^CUSTOMARILY ZEROS _____ 

C LTEMP-A HLOCK OF AT LEAST N WORDS OF TEMPORARY INTFGFR STORAGE 

C ' 

C IN ADDITION TO OVERWRITING A WITH ITS INVERSE, THE ROUTINE ALSO SETS 
C I ERR , DET , NP I V , P I V , I PR , AND LRC TO- 

C 

C 1FP.R- 0 IF INVERSION IS COMPLETED AND NO TROURLF IS OETFCTED 

C 7 IF MAGNITUDE OF CURRENT PIVOT FAILS TO FXC FED EPS 

C .... ^INVERSION WILL NOT RE COMPLFTFDS 

C 

DET-PLUS OR minus THE PRODUCT OF THE CURRENT AND ALL PREVIOUS 

PIVOTS " 

NP I V— THE NUMRFR OF THE CURRENT PIVOT *F I R ST , SFCOND, FTC . < 

P1V-THE CURRENT PIVOT 

LPR-THE FIRST NP1V POSITIONS LIST THE PIVOT ROW INDICES IN ORDER 
OF USE 

Ll’C-THF HIRST NPIV POSITIONS LIST THE PIVOT COLUMN INDICES IN 
ORDER OF USE 

DO INITIALIZATIONS " — 

1 ERR = 0 “ " ' ' ' ' 

DET= I .0 
00 7 I = 1 , N 
L PR ( I > = I 

? lpc( i ) = i " : 


7046AD 1 7 

7046A01 R 

7046A01 9 

7046A0P0 

7046A0P1 

7046A0P? 

7046A0PS 

7046A0P6 

7046A0 ? 7 

7046A0? 8 

7046A0P9 

7046A030 

7046A0 3 1 

7046A03? 

7046A033 

7046A034 

7046A03S 

7046A036 

7046A037 

7046A038 

7046A039 

7046A040 

7046A041 

7046A04? 

7046A043 

7046A044 

7046AD45 

7046AD46 

7D46A047 

7046A048 

7046A049 

7046A050 

7046A051 

7046A05? 

7046A0S3 

7046A0S4 

7046AD55 


REGIN FXCHANGF PROCESS 


00 17 NP= 1 ,N 
NP I V=N P 

SELECT PIVOT 

PIV=0.0 
DO 4 K=NP,N 

1 = L P R ( K I 

. DO 4 L-NP ,N 
J = L PC ( L I 

1FCARSI A( J,J> I-ARS(PIV) 14,3, 3 

3 K P I V = K 
L P I V = L 
I P I V= I 
JPIV-.I 

P I V = A ( I, J) 

4 CONTINUE 
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C UPDATE DETERMINATE AND PIVOT ROW AND COLUMN LISTS 
C 

DFT=DFT*P IV 
ITEMP=LPR(MP) 

LPR(NP) =LPR( KPIV) 

LPR(KPIV)=ITFMP 

]TFMP=LPC(NP) 

LPC(NP)=LPC(LPIV) 

LPC(LPIV)=ITEMP 
C ... 

C EXIT IF PIVOT TOO SMALL 

C 

I F ( E PS-ApS ( P I V ) 18,7,7 
7 I E RR = 2 
RETURN 

C 

C MODIFY PIVOT ROW 
C 

P DO 9 J=1 ,N 

9 A( IPIV,J)=-A( IPIV.JI/PIV 
C 

C _MOOIFY OTHER ROWS 

c 

DO 14 1=1, N 
IF( I-IPIV)10,14»10 
10 TEMP = A ( I , JP IV ) 

I F ( TFMP ) 1 1 , 14 , 1 1 

11 DO 13 J = 1 , N 

)F< J-JP I V ) 1 2 , 1 3 , 1 2 

12 A( I , J) =A( I , J) +A( IPIV, J )*TFMP 
1? CONTINUE 

1'- CONTINUE 
C 

C _ MODIFY PIVOT COLUMN 

C - ■ 

DO 15 I = 1 , N 

15 A ( I , JPI V ) =A( I , JP I V ) /P I V 
At IPIV, JPI V)=-A( IPIV, JP IV I 
17 CONTINUE 
C 

C END EXCHANGE PROCFSS " — 

C 

C UNSCRAMPLE ROWS OF INVERSE AMD ADJUST SIGN OF DETERMINANT 

DO IP 1=1 ,N 
L = L PR ( I ) 

IP LTEMPIL ) =LPC( I ) 

DO 22 1=1, N 

19 K = LT EMP ( I ) 

IF( I-K 120,22,20 

20 DFT = -DE T 

DO 21 J = 1 , N 

TFMP = A( I, J) — - • - 

A( I ,J)=A(K,J) 

21 A(K,»l) = TEMP 
LTEMPI I )=LTEMP(K) 

LTFMP(K)=K 

GO TO 19 

22 CONTINUE - - 

C 

c unscramble conlmns of inverse 
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DO 2? 1=1 ,M 
L = l PC ( I ) 

2? LTEMPIt. ) = LPR( 1 ) 

DO 27 1 = 1 ,N 

24 K = LT FNP (I) 

IF C 1-K 125,27, 25 

25 DO 26 0=1 ,N 
TFMP=A(J,I) 
ACJ,J)=AIJ,K) 

26 A ( J , K ) =TFMP 

11 FMP (I ) =11 F N 1 P ( K ) 
LTFMPf K ) =K 
GO TO 2 A 
27 CONTINUE 
RFTURN 
END 

/* 

//GO. SYS IN DD * 
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APPENDIX D 


ALUMINUM OXIDE X-RAY DIFFRACTION PATTERNS 
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Fig. 79 - X-ray Diffractometer Trace of Vapor Grown Aluminum Oxide 
Powder from Run VA-10 


ON 



Table XXXIII 

X-ray 

Powder Diffraction Data for Alumina, 

, Run No. VA- 

•10 



Line 

No. 

20 

d spacing 

I 

Phase 

Line 

No. 

20 

d spacing 

I 

Phase 


1 

17.2 

5.15 

M 

A0 

21 

64.4 

1.45 

M 

K, 0 


2 

19.7 

4.50 

S 

>N, K , 9 

22 

65.4 

1.43 

W 

K,0 


3 

22.2 

4.00 

M 

A 

23 

67.4 

1.39 

s 

Q!,A,7 ,N,K, 

0 

4 

25.2 

3.53 

W 

a, 6 

24 

72.6 

1.30 

w 

K 


5 

27.5 

3.24 

w 


25 

73.5 

1.29 

M 

A, 0 


6 

29.4 

3-04 

w 

K 

26 

74.3 

1.28 

M 



7 

31.5 

2.64 

s 

A,N,K ,6 

27 

75-8 

1.25 

W 

A, 0 


8 

33.1 

2.70 

M 

A,K ,0 

28 

76.9 

1.24 

W 



9 

34.6 

2.59 

S 

a,A,K,0 

29 

77.7 

1.23 

M 

N, 0 


10 

36.9 

2.43 

s 

^ ) 7 ) I* 3 K 5 0 

30 

81.9 

1.18 

W 



U 

39.1 

2.30 

M 

A,7,N,K, 0 

31 

83.9 

1.15 

w 

7, A 


12 

42.6 

2.12 

w 

a,K,0 

32 

85.6 

1.13 

M 

N 


13 

45.9 

1.9« 

M 

4j7^') ; L0 

33 

88.5 

1.10 

W 



14 

46.9 

1.94 

M 

A,K 

34 

94.1 

1.05 

M 



15 

48.2 

1.89 

M 

K, 0 

35 

96.5 

1.03 

M 

7>N 


16 

51.0 

1.79 

W 

K 

36 

99.6 

1.01 

M 



17 

57.4 

1.60 

w 

a,0 

37 

101.6 

.993 

M 



18 

60.4 

1.53 

M 

7, N,0 

38 

102.5 

.988 

W 

7 

\a5\ 

19 

61.7 

1.50 

M 

K 

39 

104.7 

.979 

M 


W|-\ 

20 

62.7 

1.48 

M 

0 

40 

108.1 

.951 

M 


\oj.O l 






41 

111.2 

.934 

M 


\<* Tl 
\ o 


I 



Fig. 80 - X-ray Diffraction Pattern for Alumina, Run VA-10 
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TITANIUM DIOXIDE X-RAY DIFFRACTION PATTERNS 
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>ING 


BL 


ANK NOT 
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Table 

XXXIV - X- 

■ray Powder Diffraction 

Data for 

Titanium, 

Run No. TiO- 

>-l 


Line 

No. 

20 

d spacing 

I 

Phase 

Line 

No. 

20 

d spacing 

I 

Phase 

1 

27.63 

3.225 

S 

Ti0 2 

19 

95.53 

1.040 

W 

Ti0 2 

2 

36.28 

2.474 

M 

Ti0 2 

20 

97.48 

1.0246 

VW 

Ti0 2 

3 

39.43 

2.283 

W 

Ti0 2 

21 

106.18 

• 9633 

VW 

Ti0 2 

4 

41.48 

2.175 

M 

Ti0 2 

22 

116.53 

.9056 

w 

Ti0 2 

5 

44.23 

2.046 

W 

Ti0 2 

23 

117.88 

.8992 

w 

Ti0 2 

6 

45.73 

1.982 

W 

Camera 

24 

120.28 

.8882 

M 

Ti0 2 

7 

54.53 

1.681 

s 

Ti0 2 

25 a 

122.88 

.8770 

M 

Ti0 2 

8 

56.78 

1.620 

M 

Ti0 2 

25b 

123.63 

.8739 

W 

Ti0 2 

9 

62.98 

1.475 

W 

Ti0 2 

26 a 

131.88 

.8459 

W 

Ti0 2 

10 

64.18 

1.450 

W 

Ti0 2 

26 b 

132.58 

.8413 

VW 

Ti0 2 

11 

69.13 

1.357 

M 

Ti0 2 

27 a 

136.68 

.8288 

W 

Ti0 2 

12 

69.93 

1.344 

W 

Ti0 2 

27 b 

137.43 

.8266 

VW 

doublet 

13 

76.63 

1.242 

vw 

Ti0 2 

28 a 

140.08 

.8195 

M 

Ti0 2 

14 

82.48 

1.168 

w 

Ti0 2 

28 b 

140.98 

.8172 

W 

doublet 

15 

84.38 

1.147 

w 

Ti0 2 

29a 

155.78 

.7879 

VW 


16 

87.78 

1.111 

vw 

Ti0 2 

30 a 

162.13 

.7797 

W 


17 

89.73 

1.092 

w 

Ti0 2 

30 b 

163.88 

.7779 

VW 


18 

90.93 

1.080 

w 

Ti0 2 









Fig. 82 - X-ray Diffraction Pattern for Titania, Run Ti0 2 -1 
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ZINC ORTHOTITANATE X-RAY DIFFRACTION PATTERNS 



15 ^ /ITT 1 


BLANtt 


173 



o 

2.55 


H 

-P~ 


Zn 2 Ti 0 4 - O 
ZnTi0 3 - A 
Ti 0 2 -Rutile* □ 


O 

1.498 


O 
1.291 

O □ 


□ 

L362 


□ 

.481 

A 


1.278 4 !34j| 1466/ 



A 

3.7 


O 

4.87 


11 O.l i **•»' 


1 

80 


60 


40 


I 

20 


- — 2d l Degrees) 


Fig. 83 - X-Ray Diffractometer Trace of Oxides from Run Zn 2 Ti0 4 -2, 
Showing Zn 2 Ti0 4 as the Major Phase 
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Table XXXV 

- X-ray 

Powder 

Diffraction Data for 

Zinc 

Orthotitanate, Run Zn 

2Ti0 4 - 

-1 

Line 





Line 





No. 

26 a 

spacing 

I 

Phase 

No. 

26 

d spacing 

1 

Phase 

1 

27.45 

3.250 

S 

Ti0 2 

24 

85.90 

1.133 

w 

Zn 2 Ti0 4 

2 

29.90 

2.989 

S 

Zn 2 Ti0 4 

25 

87.45 

1.117 

vw 

Ti0 2 ,ZnTiO_ 

3 

32.90 

2.723 

M 

ZnTi0 s 

26 

88.75 

1.103 

M 

Zn 2 Ti0 4 

4 

35.25 

2.547 

VS 

ZnTiO.^.ZnpTiO^ 

27 

89.40 

1.097 

VW 

Ti0 2 

5 

36.10 

2.489 

w 

Ti0 2 

28 

90.70 

1.084 

VW 

Ti0 2 

6 

39-20 

2.299 

vw 

Ti0 2 

29 

93.50 

1.060 

w 

Zn 2 Ti0 4 

7 

40.70 

2.218 

vw 

ZnTi0 3 

30 

95.15 

1.045 

vw 

Ti0 2 

8 

41.35 

2.185 

w 

Ti0 2 

31 

97.05 

1.030 

vw 

Ti0 2 

9 

42.80 

2.114 

w 

Zn 2 Ti0 4 

32 

101.20 

0.999 

w 

Zn 2 Ti0 4 

10 

44.10 

2.055 

vw 

Ti0 2 

33 

104.00 

0.979 

w 

Zn 2 Ti0 4 

11 

49.00 

1.860 

w 

ZnTi0 3 

34 

105.10 
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Fig. 84 - X-ray Diffraction Pattern for Zinc Orthotitanate , 
Run Zn 2 Ti0 4 -l 
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